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I.  INTRODUCTION 


In  the  early  years  of  radiation-belt  research  it  was  natural  to  focus  attention  on 
geomagnetically  trapped  protons  and  electrons,  since  these  seemed  likely  to  be  the  major 
constituents.  People  realized  that  heavier  ions  must  be  present  as  well,  but  considered 
the  presence  of  such  presumably  minor  constituents  unimportant  to  the  dynamics  of 
radiation  belts.  The  improved  satellite  instrumentation  that  has  become  available  in 
recent  years  has  drastically  altered  this  complacent  view.  Heavy  ions  are  now 
recognized  to  provide  a valuable  diagnostic  probe  for  identifying  magnetospheric  source, 
loss,  and  transport  processes.  Moreoever,  heavy  ions  are  understood  to  participate 
actively  in  naturally  occurring  magnetospheric  plasma  phenomena,  and  their  artificial 
injection  into  the  magnetosphere  in  semi-controlled  experiments  plays  a growing  role  in 
modern  space  research. 


* 

i 


This  review  is  intended  primarily  to  describe  what  has  been  learned,  and  what  can 
be  learned,  about  magnetospheric  processes  from  the  study  of  heavy  ions.  There  are,  of 
course,  additional  reasons  for  studying  heavy  ions  in  space,  e.g.,  for  understanding 
nucleosynthesis;  for  identifying  the  compositions  of  the  outer  planets  and  their  satellites 
(and  of  the  atmospheres  of  these  bodies);  and  for  analyzing  source  and  loss  budgets  for 
minor  isotopic  constituents  of  the  earth.  These  questions  are  assigned  only  secondary 
importance  here,  and  some  issues  (such  as  cosmic  rays  and  atmospheric  chemistry)  are 
not  addressed  at  all. 


Having  established  these  priorities,  we  use  the  term  "heavy  ion"  to  characterize 
any  magnetospheric  ion  other  than  a proton,  and  we  restrict  our  attention  to  kinetic 
energies  $ 10  MeV/nucleon.  The  important  heavy  ions  that  occur  naturally  in  the  earth's 


magnetosphere  are  He+,  He++,  and  0+;  ions  such  as  Ba+  and  Li+  are  important  for 
plasma-injection  experiments.  For  Jupiter  the  important  heavy  ions  are  those  of  Na  and 
S (already  observed)  as  well  as  those  of  He  and  O (and  possibly  Ne,  Mg,  and  Fe).  Ions  of 
C,  N,  and  O (in  various  charge  states)  are  important  for  both  Jupiter  and  Earth,  but  some 
instruments  now  available  cannot  distinguish  among  these  and  none  can  identify  the 
charge  state  at  energies  above  a few  keV/nucleon. 

Until  recently  magnetospheric  heavy-ion  physics  was  usually  assigned  a fairly  low 
priority  in  satellite  and  rocket  experiments,  partly  because  of  the  weight  and  complexity 
of  instrumentation  needed  to  do  a good  job  and  partly  because  of  the  popular  impression 
that  heavy  ions  had  little  to  do  with  important  magnetospheric  processes.  In  fact,  as  we 
attempt  to  show  in  this  review,  there  is  a great  deal  to  be  learned  about  auroral 
acceleration  processes,  radial  diffusion,  electrostatic  convection,  ring-current  dynamics, 
sub-auroral  red-arc  (SAR-arc)  formation,  and  magnetospheric  plasma  physics  in  general 
from  the  study  of  heavy  ions.  In  perhaps  no  other  subdiscipline  of  magnetospheric 
research  is  there  more  of  importance  to  be  learned  in  the  future,  compared  to  what  is 
already  known. 

A note  on  instrumentation  is  useful,  if  only  to  remind  us  of  two  major  gaps  in  the 
capabilities  of  instruments  that  have  so  far  been  flown.  There  are  four  basic  types  of 
instrument,  some  of  which  can  be  combined:  (1)  threshold  detectors  (TD's),  usually 
made  of  silicon  (Si),  which  discriminate  heavy  ions  by  their  greater  dE/dx  (energy 
deposition  per  unit  path  length)  and  which  have  an  effective  threshold  of  a few  hundred 
keV/nucleon  for  heavy  ions;  (2)  E- dE/dx  detectors,  which  measure  dE/dx  in  a thin 
detector  and  E in  a thick  detector  which  stops  the  particle,  having  a threshold  of  a few 
hundred  keV/nucleon;  (3)  electrostatic  analyzers  (ESA's),  followed  by  momentum  or 
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velocity  filters,  which  have  an  effective  upper  cutoff  - 20  keV/charge;  and  (4)  time-of- 
flight  detectors  (TOFD's),  which  must  be  combined  with  instruments  of  other  types  in 
order  to  identify  heavy  ions.  No  data  have  yet  been  pun  ished  from  the  few  operational 
satellite  experiments  that  use  TOFD's.  The  other  instrument  combinations  that  have 
been  flown  cannot  be  used  in  the  range  from  20  keV/charge  to  100  keV/nucleon 
(approximately),  and  the  identification  of  charge  states  has  been  contingent  on 
restrictive  assumptions  about  the  charge  and  mass  composition  of  the  ion  sample 
observed.  The  advantage  of  threshold  detectors  (TD's)  is  their  simplicity.  This  accounts 
for  their  common  use  in  space.  Being  one-parameter  identification  systems,  however, 
TD's  have  less  mass  resolution  and  a greater  susceptibility  to  background  pile-up  than  the 
other  (more  complicated)  detector  systems.  Thus,  TD's  cannot  reliably  distinguish  among 
C,  N,  and  O;  moreover,  because  of  background  effects,  they  are  not  useful  if  the  flux  of 
heavy  ions  is  less  than  ~10-4  that  of  protons. 

Despite  some  of  these  observational  constraints,  significant  progress  has  been 
made  in  the  study  of  magnetospheric  heavy  ions,  and  the  future  for  research  in  this  area 
is  a very  bright  one  indeed.  We  have  organized  the  remainder  of  this  review  into  sections 
that  (in  order)  deal  with  auroral-magnetospheric  processes;  heavy-ion  tracers  for  radial 
diffusion  and  electrostatic  convection;  multi-ion  plasma  physics;  artificial  plasma- 
injection  experiments;  and  heavy-ion  physics  at  Jupiter.  Finally,  a program  for  the 
future  is  outlined  in  Section  VII. 
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II.  AURORAL-MAGNETOSPHERIC  COUPLING  AND  THE  0+  PROBLEM 

141 

In  1972  Shelley  et  al.  reported  a truly  shocking  discovery  based  on  their  mea- 
surements from  a low-altitude  polar-orbiting  satellite,  viz.,  large  fluxes  of  precipitating 

O*  ions  having  energies  up  to  12  keV  (the  upper  limit  of  their  instrument).  These  fluxes, 

2 

ranging  up  to  0.4  erg/cm  -sec-ster  in  energy  content,  were  observed  on  the  night  side  at 

L values  as  low  as  2.5  and  sometimes  exceeded  proton  precipitation  fluxes  at  the  same 

energy.  Since  particles  having  A/Z  (mass-to-charge  ratio)  = 16  in  the  solar  wind  are  only 
-3  8 

~ 10  as  abundant  as  protons  there  , it  is  fairly  clear  that  the  observed  heavy  ions  were 
of  ionospheric  origin. 

Further  studies*'*®  139,142,143  ^ave  reveaie(j  ^at  the  observed  0+  ions  generally 
precipitate  equatorward  of  the  precipitating  protons  that  accompany  them  (see  Figure  1). 
The  invariant  latitudes  of  the  0+  and  H+  precipitation  peaks,  i.e.,  their  values  of 
cos“*(L-^),  are  well  correlated  with  each  other  and  with  the  invariant  latitude  of  the 
plasmapause,  which  lies  equatorward  of  both  precipitation  maxima.  Moreover,  the 
precipitating  0+  intensity  is  well  correlated  in  time  with  the  ring-current  intensity  (as 
measured  by  the  geomagnetic  index  Dgt)  and  with  substorm  activity  (as  measured  by  the 
index  AE).  The  possibility  thus  exists  that  0+  is  a significant  contributor  to  the  storm- 
time ,-:ng  current.  Unfortunately  there  are  as  yet  few  trapped-particle  measurements 
that  distinguish  among  ion  species  in  the  ring  current. 

In  principle,  however,  it  is  possible  to  decide  whether  there  is  indeed  a significant 
0+  flux  among  the  ring  current  ions.  Because  of  the  fourfold  difference  in  number 
density  required  to  produce  the  same  flux,  a given  flux  of  0+  ions  at  any  given  energy 
would  contribute  four  times  as  much  to  Dgt  as  would  the  same  flux  of  protons;  the  0+ 
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Figure  1.  Precipitating  energy  fluxes,  nightside  (~  0300  LT)  and  dayside  (~1500  LT), 

137 

associated  with  two  magnetic  storms 
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and  H+  ions  would  differ  by  a factor  of  four  in  velocity,  but  their  drift  rates  would  be  the 

same.  Thus,  the  values  calculated  from  ion  fluxes  observed  at  the  equatorial 

satellite  S (Explorer  45)  under  the  assumption  that  these  ions  were  protons  would  yield 

discrepancies  with  Dgt  values,  measured  at  the  network  of  ground  stations  during  the 

same  time  period,  if  there  were  a significant  contribution  from  0+  in  the  ring-current  ion 

flux.  For  example,  a 5%  contribution  from  0+  in  the  observed  ion  flux  would  result  in  a 

15%  greater  ion  Dgt  on  the  ground  than  one  would  calculate  under  the  assumption  of  a 

100%-proton  flux.  Of  course,  one  cannot  neglect  the  contribution  of  ring-current 

12 

electrons  to  at  this  level.  Berko  et  al.  have  carried  out  such  a comparison  and 
have  concluded  that,  during  a time  period  when  the  Lockheed  ESA  detected  significant 
precipitation  of  0+,  at  most  a few  percent  of  the  ring-current  ion  flux  could  have  arisen 
from  0+.  In  view  of  the  compounded  uncertainties  inherent  in  such  a comparison, 
however,  it  would  be  wise  to  wait  for  direct  ion-composition  measurements  in  the  ring 
current  before  closing  the  door  on  this  question. 

The  discovery  of  energetic  0+  was  so  shocking  for  two  reasons.  This  was  the  first 

reported  example  of  a heavy  ion  as  the  dominant  component  of  a magnetospheric  plasma, 

and  no  one  at  the  time  could  think  of  a good  mechanism  for  energizing  such  heavy  ions. 

The  idea  of  cyclotron-resonant  energy  transfer  (see  Section  IV)  from  protons  to  0+  has 
22  30 

been  proposed  ’ but  is  not  very  promising.  The  idea  of  charge  exchange  between 
energetic  protons  and  atomic  oxygen  is  not  very  appealing  either,  since  this  would  tend 
to  produce  energetic  neutral  H and  cold  0+.  At  the  time  the  energetic  0+  was 
discovered,  there  had  been  no  direct  observations  of  the  parallel  (to  B)  electric  fields 
that  might  accelerate  thermal  ionospheric  0+  to  energies  of  several  keV.  Such  parallel 
C fields  had  been  contemplated  theoretically  • * and  were  invoked  from  time  to 

MM 
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time  in  the  literature,  but  the  truly  compelling  argument  for  their  existence  was 
provided  by  Evans'*1  in  his  elegant  analysis  of  precipitating  electrons  and  their 
backscatter  spectrum. 

Moreover,  the  observational  situation  has  greatly  improved  since  the  launch  of  the 
S3-3  satellite  into  an  elliptical  orbit  with  polar-region  apogee  («8000  km)  and  perigee 
( - 240  km).  Using  data  from  their  triaxial  £-field  antenna  on  S3-3,  Mozer  etal.115  have 
reported  perpendicular  (to  B)  electric  fields  as  strong  as  400  V/km  at  an  altitude 

•MIK 

- 7000  km  and  parallel  (to  B)  electric  fields  with  the  smaller,  but  stiil  unbelievably  high, 

MM 

strength  of  100-200  V/km.  An  example  of  their  observational  results,  after  spin- 

demodulation  of  the  £-field  data,  is  shown  in  Figure  2.  The  perpendicular  £ field 

showed  a distinctive  sign  reversal,  such  as  would  be  produced  by  a sheet  of  negative 

charge,  over  a latitudinal  extent  '-100  km  at  the  sites  of  enhancement  in  £n  (see 

112  144 

Figure  2).  This  feature  is  seen  ’ in  association  with  beams  of  downgoing  electrons 
and  upgoing  ions  (as  well  as  with  unbeamed  particles)  having  energies  that  range  from 
several  hundred  eV  to  a few  keV,  as  is  illustrated  by  the  gray-scaled  energy-time 
spectrogram  shown  in  Figure  3 (provided  by  P.  F.  Mizera  from  published11**  and 
unpublished  material,  1977).  This  figure  shows  an  "inverted-V"  structure  in  the  electron 
data  (i.e.,  a rise  and  fall  in  the  mean  electron  energy)  quite  similar  to  to  those  reported 

cq  en 

by  Frank  and  Ackerson  * . Electron  beams  are  not  obvious  here,  because  of 

decollimation  by  the  earth's  magnetic  field  and  contamination  by  backscatter.  However, 
the  ion  beams  are  very  distinct,  partly  because  the  earth's  ^ field  collimates  them.  A 
measure  of  the  ionic  composition  of  one  such  beam  is  shown  in  Figure  4.  The  virtual 
absence  of  upgoing  helium  ions  is  evident  here.  Further  analysis  reveals  that  the  H+  and 

+ 8 —9  — 1 .1  _ 1 

O beams  are  of  comparable  intensity,  ranging  up  to  ~ 10  cm  sec  ster  keV  at  the 
144 

peak,  which  occurs  at  an  energy  ~ 1 keV. 
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Figure  2.  Electric-field  components  observed  at  S3-3  on  poleward-bound  pass  through 
northern  auroral  zone11'*,  UT  = 41.08—41.30  ksec. 
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Such  results  are  readily  interpreted  in  terms  of  upward  parallel  (to  B)  electric 

MW 

112 

fields;  it  appears  in  at  least  one  case  that  the  potential  drop  below  the  7000-km 

satellite  altitude  totalled  ' 2 kV  and  that  the  potential  drop  above  the  satellite  totalled 

an  additional  1 keV  or  so.  The  corresponding  parallel  £ fields  in  this  case  would  amount 

115 

to  - 0.3  V/km,  or  about  0.2%  of  those  reported  by  Mozer  el  al.  . Theoretical  expecta- 

99  70  151  152 

tions  would  place  e < 1 V/km.  Various  proposed  configurations  ’ ’ for  the 

equipotential  contours  in  such  an  auroral  feature  are  illustrated  in  Figure  5.  For 

example,  Swift  has  suggested  a nearly  parallel  (to  ^B)  electrostatic  "shock" 

configuration  with  a spatial  scale  of  several  ion  Larmor  radii  (tens  of  km)  in  latitude  (see 

Figure  5c).  Such  a configuration  could  account  for  the  observed  reversal  in  £^  and  for 

"inverted-V"  particle  signatures.  In  Swift's  model*^  the  "shock"  is  driven  by  energetic 

ions  (other  than  those  in  the  beams)  and  is  essentially  an  electrostatic  ion-cyclotron  wave 

(see  Section  IV)  at  zero  frequency.  Additional  sources  of  parallel  electric  field  include 

1 29 

electrons  and  ions  with  different  pitch-angle  anisotropies  ’ and  anomalous  resistivity, 
which  may  furnish  the  dissipation  needed  to  convert  Swift's  structure  into  a true  shock. 

What  happens  to  the  upgoing  ions  after  they  pass  through  the  S3-3  altitude?  There 

are  few,  if  any  downgoing  ion  beams,  as  conjugacy  might  allow,  but  similar  ions  beams 

have  been  seen  at  the  equator1  , and  streaming  O has  been  observed  in  the  earth's 

6 1 

magnetotail  at  energies  of  1-5  keV.  The  latter  observation  would  seem  to  require  an 

8 -2  -1  112 
ionospheric  source  flux  -10  cm  sec  , which  is  comparable  to  that  observed  at 

S3-3.  Perhaps  conjugate  upgoing  beams  thermalize  each  other  through  a two-stream 

133 

instability  as  they  reach  the  equator  in  closed  field  lines  while  upgoing  streams  on 

9 

open  field  lines  escape  from  the  magnetosphere  in  the  manner  of  polar-wind  ions  . 


Reconnection  and  convection  may  introduce  some  ambiguity  in  the  above  alternatives. 
Ions  of  energy  1-5  keV  that  become  trapped  on  closed  field  lines  in  consequence  of  such 


processes  would  become  part  of  the  ring  current.  Those  with  mirror  points  near  the 

-3 


equator  are  subject  to  energization  in  proportion  to  L through  sunward  convection  or 
radial  diffusion.  Thus,  ions  of  energy  10-50  keV  might  be  seen  at  L values  near  the 


plasmapause  (but  not  far  inside,  since  the  plasmapause  is  an  approximate  boundary  for 


low-energy117  and  medium-energy40  plasma  access).  Oxygen  ions  in  the  ring  current 


would  charge-exchange  at  only  ~20%  of  the  proton  rate1  0 for  E<10  keV,  and  so  their 
fraction  of  the  ring  current  at  low  energies  would  increase  with  time  until  a steady  state 
is  attained  (see  Section  IV). 


The  foregoing  scenario  represents  one  school  of  thought  concerning  the  origin  of 


ring-current,  radiation-belt,  and  auroral  ions.  In  the  sense  that  it  points  toward  an 

, o 

ionospheric  source,  this  scenario  resembles  the  proposal  that  the  polar  wind  provides  the 

plasma  sheet  with  its  charged  particles.  Of  course,  the  electric  field  associated  with  the 

polar  wind  is  much  weaker  than  that  associated  with  an  "inverted-V"  structure,  but  the 

latitudinal  extent  of  the  polar  wind  is  much  greater.  A major  difference  related  to 

electric-field  strength,  however,  is  that  only  the  "inverted-V"  structure  can  supply 

significant  0+  to  the  plasma  sheet. 

The  other  school  of  thought  concerning  the  origin  of  magnetospheric  ions  is  that 

58 

they  enter  through  the  dayside  cleft  and  through  the  flanks  of  the  nightside  magneto- 
pause. The  dayside  cleft  is  associated  with  the  regions  of  weak  B near  the  neutral 
points100.  The  nightside  magnetosphere  is  separated  into  two  lobes  by  the  equatorial 
current  sheet,  which  occupies  a region  of  weak  magnetic  field  comparable  to  that  of  the 
dayside  cleft.  Charged  particles  that  enter  the  magnetosphere  from  the  outside  are 
mainly  of  solar  origin. 


6 

The  above  considerations  represent  a major  dilemma  in  magnetospheric  physics: 
Do  energetic  particles  in  the  ring  current,  radiation  belts,  and  aurora  have  their  main 

g 

source  in  the  ionosphere  or  in  the  solar  wind?  Axford  had  proposed  isotopic  and  charge- 

state  observations  on  helium  ions  in  the  magnetosphere  as  a way  to  identify  the  dominant 

source,  since  the  He++/He+  and  the  He^/He^  abundance  ratios  are  much  greater  in  the 

solar  wind  than  in  the  ionosphere.  Such  measurements  have  since  been  made  in  diffuse 

proton  auroras  (which  are  quite  different  from  the  well  localized  visual  auroral  arcs 

associated  with  "inverted-V"  structures  in  electron  spectrograms  and  with  upgoing  ion 

beams)  and  also  during  auroral  break-up.  A 1973  review  of  such  measurements  has  been 

provided  by  Reasoner  , who  (following  Cornwall  0 ) has  indicated  the  need  for 

caution  in  interpreting  charge-state  results  when  atmospheric  collisions  may  have 

170  171 

intervened  between  the  source  and  the  observer  (see  Section  III).  Whalen  etal.  ’ 

used  an  ESA  (see  Section  I)  and  momentum  filter  (magnetic  deflection  system)  on  a 

rocket  flight  to  distinguish  between  He+  and  He++  among  precipitating  ions  in  diffuse 

proton  auroras.  They  found  no  statistically  significant  abundance  of  He+  compared  to 

He++  and  therefore  concluded  that  the  source  was  solar  rather  than  ionospheric.  The 

5 -9  -1  -1 

precipitating  ion  flux  was  ~ 10  cm  sec  ster  keV,  or  less  than  ~1%  of  the  upgoing  ion 

flux  found  in  association  with  "inverted-V"  structures.  Similar  results  were  found  by 
135 

Sharp  et  al.  in  satellite  passes  through  the  auroral  zone.  Recoverable  metal  foils 

have  been  flown  on  auroral  rockets  and  analyzed  with  a mass  spectrometer  to  ascertain 
3 4 7 

the  He  /He  abundance  ratio  , which  turned  out  during  an  auroral  break-up  to  approxi- 

27 

mate  that  of  the  solar  wind  . These  results  provided  additional  support  for  the 
hypothesis  of  a solar  (rather  than  ionospheric)  source. 
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However,  He  has  been  found  precipitating  in  the  auroral  zone  during  a storm  at 
7 -2  -1  -1  -1 

a flux  level  ~ 10  cm  sec  ster  keV  for  energies  of  a few  keV.  This  level  exceeded 

that  of  the  precipitating  H+  and  0+  fluxes  in  an  auroral  event  that  yielded  no  detectable 

He++.  The  ions  seen  in  this  instance  were  probably  not  of  solar-wind  origin:  both  the 

He+/He++  and  the  0+/H+  ratios  were  characteristic  of  the  ionosphere.  Thus,  the 

observational  situation  seems  to  be  somewhat  variable.  Further  opportunities  to  identify 

the  source  of  auroral,  ring-current,  and  radiation-belt  ions  are  potentially  available  from 

a study  of  ions  heavier  than  helium.  Krimigis  etal.  , for  example,  have  found  that  ions 
-3 

having  Z?  3 are  ~10  as  abundant  as  helium  ions  having  the  same  energy/nucleon.  The 
13 

challenge  is  to  distinguish  such  heavy  ions  (A  i 12)  from  each  other  at  E/A  S 1 
MeV/nucleon;  determination  of  the  C/O  abundance  ratio  in  this  range  should  provide  the 

most  decisive  possible  test  for  identifying  the  particle  source,  since  this  ratio  is  larger  by 

5 13 

at  least  a factor  ~ 10  in  the  solar  wind  than  in  the  ionosphere  . 

1 

. 
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III.  HEAVY  IONS  AS  PROBES  FOR  MAGNETOSPHERIC  PROCESSES 


Before  the  discovery  of  energetic  0+  ions  in  great  numbers141,  the  only  known 
role  for  heavy  ions  in  magnetospheric  processes  was  to  serve  as  test  particles  to  probe 
the  dynamics  of  protons.  By  studying  heavy  ions  one  hoped  to  discover  the  primary 
source  of  radiation-belt  particles,  as  well  as  the  mechanisms  by  which  they  are  injected 
into  the  region  of  closed  field  lines,  transported  across  L,  and  energized.  Of  course, 
heavy  ions  are  now  of  intrinsic  interest  also,  aside  from  their  historic  role  as  test 
particles. 

The  earliest  quantitative  observations  of  trapped  helium  ions  are  those  reviewed 
92  93  1 fi 

by  Krimigis  ’ and  by  Blake  and  Paulikas  . The  data  were  taken  from  low-altitude 

satellites  in  polar  orbits,  i.e.,  from  a TD  on  Injun  5 and  from  an  E -dE/dx  telescope  on  the 

OV1-19  satellite.  The  orbits  precluded  the  generation  of  a profile  in  L for  particles 

1 7 57  94 

mirroring  near  the  equator.  Many  later  observations  shared  this  drawback  ’ ’ . 
Representative  off-equatorial  proton  and  helium-ion  spectra  are  shown  in  Figure  6 
with  the  indicated  normalization.  The  instrument  was  unable  to  distinguish  He+  from 
He++;  the  term  "alpha  particles",  as  commonly  used  in  this  context,  includes  both  charge 
states.  The  proton  and  alpha  spectra  corresponding  to  the  same  L and  B/Bq  (ratio  of 
local  to  equatorial  field  intensity  on  a given  field  line)  in  Figure  6 seem  to  obey  the  same 
asymptotic  power  law.  The  index  in  the  power  law  would  therefore  be  unaltered  if  one 
plotted  flux  against  energy  per  nucleon  (E/A)  or  energy  per  charge  (E/Z),  rather  than 
energy  per  entity  (E)  as  is  done  in  Figure  6.  However,  the  normalization  would  be 
altered.  A comparison  of  these  three  conventions  as  they  affect  the  ratio  J^/J^  is 
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Figure  6.  Differential  alpha-particle  and  renormalized  proton  spectra  from  OV1-14 
(1968-26B)  at  selected  off-equatorial  locations  (from  B/Bq  = 2.28  at  L = 2.3  to  B/BQ  = 
8.48  at  L = 3.3).  The  indicated  renormalization  factors  represent  the  o/p  ratio  at  fixed 
E,  rather  than  at  fixed  E/A. 
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shown  in  Figure  7 for  alpha-particle  energies  in  the  range  1-8  MeV.  A certain 

159 

dynamical  significance  corresponds  to  each  convention.  For  example,  Tverskoi  and 
79 

Hess  had  predicted,  on  the  basis  of  a theory  to  be  discussed  below,  that  the  a/p  flux 

81 

ratio  at  fixed  energy/nucleon  (E/A)  should  be  near  the  solar-wind  ratio,  which  is  found 
to  be  ~ 4%.  The  actual  a/p  ratio  at  fixed  E/A  turns  out  to  be  -100  times  smaller  than 
this,  at  least  for  the  off-equatorial  observations  summarized  in  Figure  7. 

qc  -l  co 

Analysis  of  quiet-time  Injun-5  measurements  on  CNO  ions  ’ revealed  a 

_3 

similar  puzzle:  the  CNO/He  ratio  at  E/A  > 0.3  MeV/nucleon  was  only  ~ 10  in  the 

_o 

magnetosphere,  as  compared  with  ~10  in  either  the  solar  wind  or  the  ionosphere. 

Thus,  the  CNO/H  ratio  is  -1000  times  smaller  in  the  off-equatorial  magnetosphere 

13 

than  in  the  solar  wind  during  quiet  times.  It  would  have  been  extremely  valuable  to 

distinguish  among  the  separate  species  C,  N,  and  O,  but  this  was  not  possible  with  the  TD 

aboard  Injun  5.  Carbon  and  oxygen  have  been  measured  separately  in  the  radiation  belts 

at  much  higher  energies  (13-33  MeV/nucleon).  The  resulting  quiet-time  O/C  ratio  of 

0.5  ± 0.4  is  consistent  only  with  an  extraterrestrial  source  . Unfortunately,  it  has  been 

163 

impossible  to  establish  the  CNO/He  ratio  reliably.  Verzariu  has  reported  observing  a 
_2 

CNO/He  ratio  -10  in  the  magnetosphere  following  a solar-particle  event  in  which  the 
magnetospheric  CNO  flux  increased  by  a factor  - 100  compared  to  quiet-time  levels. 

The  general  impression  of  having  too  few  magnetospheric  heavy  ions  was  changed 
considerably  when  equatorial  measurements  became  available1  It  then  became 

clear  that  the  equatorial  pitch-angle  distribution  of  helium  ions  is  considerably  narrower 
than  that  of  protons,  which  would  imply  that  the  a/  p ratio  is  a sharply  decreasing 
function  of  magnetic  latitude.  An  example1'*  from  the  inner  zone  (Ls2)  is  shown  in 
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Figure  7.  Observed  a/ p ratios  at  fixed  E,  fixed  E/Z,  and  fixed  E/ \.  magnetic 
latitude  of  observation  increases  strongly  with  L. 


-26- 


I 


Figure  8,  where  a fit  to  distributions  of  the  form  sinnaQ  yielded  n = 6.7  for  protons  and 

64 

n = 11.3  for  helium  ions.  Fritz  and  Williams  report  ns:  8 for  outer-zone  (3  s L <5)  alpha 

particles  and  n « 5 for  outer-zone  protons  at  E/A  « 0.2-0. 5 MeV/nucleon.  They  find 

-2 

o/p  ratios  approaching  10  at  the  magnetic  equator,  in  rough  agreement  with  solar- 
wind  and  ionospheric  values.  Figure  9 shows  that  the  pronounced  variation  of  the  at p 
flux  ratio  with  equatorial  pitch  angle  (hence,  magnetic  latitude)  is  rather  insensitive  to  L 

56 

over  a certain  range  of  L values  . It  would  have  been  difficult  to  sort  out  any  such 
variation  from  Figure  7 (above),  since  B/Bq  there  was  an  increasing  function  of  L rather 
than  a constant  parameter  of  the  profile. 

There  exist  major  temporal  variations  in  the  a/p  ratio.  Magnetic  storms  and 
substorms  cause  it  to  increase  by  as  much  as  an  order  of  magnitude  far  off  the  equator, 
but  by  far  less  at  the  equator  itself.  For  E/A  > 0.2  MeV/nucleon  the  alpha-particle 
intensity  far  off  the  equator  decays  more  slowly  than  would  be  predicted  from  the  theory 
of  ionic  collisions  with  the  atmosphere  (which  theory  entails  primarily  energy  deposition 
and  charge  exchange)  following  a magnetic  storm.  This  finding  suggests  a role  for  pitch- 
angle  diffusion  by  wave-particle  interactions  (see  Section  IV).  It  has  already  been 
mentioned  that  magnetospheric  CNO  fluxes  may  increase  by  two  orders  of  magnitude 
during  a magnetic  storm.  Evidently  the  injection,  transport,  energization,  and  loss 
process  depend  in  important  ways  on  Z and  A. 

79  159 

Until  1971  the  only  theoretical  framework  available  ’ for  analyzing  such 
results  for  E/A  2 0.2  MeV/nucleon  was  one  that  postulated  a universal  radial  diffusion 
coefficient  (independent  of  E,  Z,  and  A)  and  a sole  source  originating  in  the  solar 
wind.  The  radial  diffusion  was  supposed  to  be  caused  by  magnetic  impulses  that  are  step- 
like on  the  drift  time  scale.  It  was  assumed  that  injection  into  the  magnetosphere  at 
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Figure  8.  Equatorial  pitch-angle  distributions  of  alpha  particles  (left  panel)  and 
protons  (right  panel)  from  OV1-19  (1969-26C)  data1!\ 


Figure  9.  Observed  a/p  flux  ratios  at  fixed  E/A,  equivalent  to  J^/J^  at  B/B 
esc2  a , as  compiled56  from  the  data  of  various  spacecraft  instruments. 


L2  7 somehow  separately  "thermalized"  the  various  solar  wind  ions  so  as  to  produce  a 

temperature  proportional  to  A.  In  this  case  the  source  spectrum  would  have  the  same 

He++/H+  ratio  at  each  E/A,  viz.,  the  a/ p abundance  ratio  in  the  solar  wind.  Moreover, 

the  universality  of  the  radial-diffusion  coefficient  would  preserve  this  ratio 

throughout  the  magnetosphere  despite  the  presence  of  Coulomb  interactions  (which 

o 

deposit  E/A  at  a rate  solely  proportional  to  Z /A  for  ions  having  the  same  velocity). 

79  159 

Charge  exchange  is  neglected  ’ at  the  energies  of  interest  but  would  likely  erode  the 

proton  distribution.  The  off-equatorial  measurements  of  the  a/p  flux  ratio  are  in 

substantial  disagreement  with  the  major  prediction  of  this  theory,  viz.,  that 
-2 

>10  throughout  the  magnetosphere. 

The  above  discrepancy  can  be  alleviated  by  postulating  a partial  or  complete 

mutual  "thermalization"  of  solar-wind  ions  to  produce  the  source  spectra.  In  this  case 

one  might  expect  the  alpha-particle  source  spectrum  to  have  a temperature  less  than 

four  times  that  of  the  protons.  This  is  a reasonable  model  if  the  protons  and  heavy  ions 

enter  the  trapping  region  via  the  turbulent  nightside  plasma  sheet  rather  than  directly 

from  the  bow  shock.  While  it  is  true  that  all  solar  wind  ions  have  about  the  same 
81 

velocity  , and  hence  a mean  energy  proportional  to  A,  their  entry  via  the  plasma  sheet 

could  cause  them  to  share  the  available  energy  so  as  to  produce  a mutual 
3C 

thermalization  In  this  case  one  should  expect  the  a/p  ratio  of  the  source  spectra  to 

be  an  exponentially  decreasing  function  of  E/A,  which,  moreover,  would  be  mapped 

1 S9 

throughout  the  equatorial  magnetosphere  in  accordance  with  the  theory  of  Tverskoi 
79 

and  Hess  . Mutual  thermalization  would  also  account  for  the  observational  fact  that 
alpha-particle  and  proton  pitch-angle  distributions  are  differently  anisotropic  in  the  inner 
magnetosphere  (see  Figures  8-9,  above).  The  source  of  extreme  anisotropy  in  either 


species  resides  in  the  fact  that  radial  diffusion  from  an  external  source  energizes 
particles  by  a factor  that  varies  inversely  with  mirror  latitude1  ^ Thus,  as  compared 
with  particles  that  mirror  on  the  equator,  those  that  mirror  at  higher  latitude  must  have 
enjoyed  a higher  E/A  in  the  source  spectrum  in  order  to  have  the  same  energy /nucleon  at 
an  interior  L value.  Thus,  the  steepness  of  the  source  spectrum  manifests  itself  as  an 
anisotropy  in  the  equatorial  pitch-angle  distribution  at  interior  L values.  This  is  true  for 
either  species.  However,  with  partial  or  complete  mutual  "thermalization"  accom- 
plished, the  source  spectrum  for  alpha  particles  is  steeper  with  respect  to  E/A  than  the 

source  spectrum  for  protons.  Therefore,  the  alpha  particles  show  an  even  larger 

131 

anisotropy  at  interior  L values  than  do  the  protons 

The  actual  situation  is  much  more  complicated  than  this.  Ions  having  E/A  ~ 0.2 
MeV/nucleon  at  L = 2-3  may  have  had  energies  as  low  as  20  keV/nucleon  in  the  source 
spectrum  at  the  site  of  injection  (nominally  L = 7);  if  the  first  two  adiabatic  invariants 

_3 

are  conserved  in  the  process,  as  is  assumed,  then  the  nonrelativistic  energy  varies  as  L 

for  a - 90°  and  as  L 2 in  the  unattainable  limit  o = 0.  The  variation  of  E/A  with  L is 
o o 

more  complicated  at  intermediate  values  of  oq,  since  oQ  itself  then  varies  with  L.  At 
energies  such  as  20  keV/nucleon  it  turns  out  that  charge  exchange  is  a very  important 
loss  (neutralization)  process  that  depends  strongly  on  Z and  A;  moreover,  charge 
exchange  efficiently  converts  He  to  He  at  such  low  energies  ’ , as  can  be  seen 

from  Figure  10.  One  also  expects  wave-particle  interactions  (see  Section  IV)  to  cause 
pitch-angle  diffusion  of  both  protons  and  alpha  particles,  but  wave-particle  interactions 
can  act  with  vastly  different  strengths  on  the  various  ionic  species  if  the  interactions  are 
cyclotron-resonant.  Finally,  radial  diffusion  by  fluctuations  of  the  electric  "convection” 
field,  such  as  might  occur  during  substorms,  is  a very  important  process  to  consider,  and 
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its  coefficient  has  a strong  dependence  on  the  ratio  of  E/A  to  Z/A,i.e.,  on  E/Z.  It  was 

o J 

shown  by  Cornwall  that  such  electrostatic  radial  diffusion  produced  a somewhat  better 

fit  to  the  proton  data  then  available  than  did  the  type  of  radial  diffusion  caused  by  step- 

69 

like  magnetic  impulses,  and  by  Gregory  that  charge  exchange  could  strongly  affect 

35 

pitch-angle  distributions.  Cornwall  subsequently  attempted  to  include  most  of  the 
above  effects  (as  well  as  Coulomb  loss  and  magnetic  radial  diffusion)  in  a detailed 
calculation  of  the  magnetospheric  transport  and  loss  of  helium  ions  and  protons. 

The  electrostatic  and  magnetic  diffusion  coefficients  are  each  proportional  to  a 

spectral  density  of  c£/B,  evaluated  at  the  azimuthal-drift  frequency  of  the  particle  in 
54 

question  . The  aziruthally  asymmetric  component  of  the  £ field  induced  bv  step-like 
changes  in  J3^has  a flat  spectrum,  since  the  corresponding  3J3/at  consists  of  a random 
series  of  Dirac  delta  functions  in  time.  Moreoever,  the  azimuthally  asymmetric 
component  of  this  £ varies  as  r , since  its  curl  must  be  proportional  to  time  derivative 

MM 

of  a B-field  perturbation  that  increases  linearly105  with  the  geocentric  distance  r.  Since 

O 

the  main  (dipolar)  J3  field  in  c£/B  varies  as  r , the  result  is  a radial-diffusion 
coefficient  proportional  to  L10  and  independent  of  drift  frequency  (since  the 

spectrum  of  £ is  flat).  A more  detailed  analysis 54,91  designed  to  account  for  the 

/ \ 

correct  mapping  of  £ and  B throughout  the  magnetosphere  introduces  in  D; , a factor 

MM  AMP  L Ll 

that  depends  strongly  upon  the  mirror  latitude  of  the  particle,  i.e.,  upon  aQ.  The 

(e) 

electrostatic  impulses  that  generate  D'L£  are  assumed  to  result  from  an  equatorially 
uniform  £ field  that  rises  sharply  at  random  times  and  subsequently  relaxes  toward  its 

MM 

mean  value  with  a decay  time  t ~20  min.  The  corresponding  spectrum  is  flat  for  drift 

2 

periods  2ir/w^»2hr  and  inversely  proportional  to  w ^ for  drift  periods 

o 

2ir/u>D  » 2 hr.  The  drift  frequency  itself  is  proportional  to  M/L  Z for  particles  of 
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first  adiabatic  invariant  (magnetic  moment)  M that  mirror  near  the  magnetic  equator, 

and  there  is  a weakly  inverse  variation  of  with  mirror  latitude  for  particles  having 

the  same  kinetic  energy  E.  Thus,  the  spectral  density  of  c£/B  at  is  roughly 

proportional  to  L6  for  2Wu>D  » 2 hr  and  to  L19  (Z/M)2  sin4^  for  2tt  / <*' ^ « 2 hr. 
35 

Cornwall  has  summarized  these  considerations  by  setting 


(1) 

for  particles  that  mirror  near  the  equator,  where  is  the  kinetic  energy  (in  keV)  at 
L = 7 (or  equivalently,  the  magnetic  moment  in  MeV/G)  and  1.3  s KS13.  These  values 
have  subsequently  been  endorsed  (with  K = 2.6)  by  Spjeldvik148  in  his  analysis  of  proton 
data,  but  other  investigators22,48,85  prefer  a much  larger  magnitude  for  D^,  perhaps 
5-10  x 10-9  L10  day'1. 

For  particles  having  the  same  value  of  E?/A  ( > 40  keV/nucleon)  it  is  evident  from 
(1)  that  is  roughly  proportional  to  L10(Z/A)2  for  L S 7.  Thus,  the  electrostatic 

radial  diffusion  coefficients  for  H+:He+*:He+  bear  the  ratios  16:4:1.  Since  charge 
exchange  acts  quite  rapidly  (on  the  diffusion  time  scale)  to  convert  He++  into  He+  (see 
Section  IV),  the  net  result  is  that  helium  ions  diffuse  much  more  slowly  from  the  external 
(L  * 7)  source  than  do  protons,  and  so  the  helium  ions  have  (on  the  average)  more  time  to 
be  lost  (from  whatever  cause)  than  do  protons  before  reaching  the  inner  magnetosphere 
(L  * 2-3).  This  condition  is  somewhat  mitigated  by  the  fact  that  loss  rates  due  to  charge 


} 


D = 

LL  ULL  u LL 


= ( 10'5  K [L4  + (t/20  min)2  (E?/Z)2  j _1 


+ 2 x 10"10  } L10  day"1 
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exchange  (He+  — He°)  and  Coulomb  interactions  (deposition  rate  of  E/A  proportonal  to 

2 + 

Z /A)  are  larger  for  protons  than  for  He  at  the  same  E/A. 

35 

Figures  11  and  12  show  some  of  the  result  of  Cornwall's  numerical  calculation  . 

Figure  11  illustrates  how  efficiently  charge  exchange  has  increased  the  He+/He++  ratio 

from  zero  at  the  source  (L  = 7)  to  ~1  at  L < 4.  The  calculation  has  included  the  cross 

section  for  He+— He++,  which  overtakes  the  other  two  helium-ion  charge-exchange  cross 
31 

sections  at  a kinetic  energy  ~1  MeV.  In  order  to  avoid  economic  difficulties  that 

35 

seemed  to  threaten  at  the  time,  Cornwall  considered  only  particles  of  vanishing  second 

invariant  J (i.e.,  particles  mirroring  at  the  magnetic  equator)  in  his  numerical  work.  This 

was  unfortunate,  since  the  only  observed  «/p  ratios  then  available  for  comparison  had 

-4 

been  obtained  far  off  the  equator.  These  were  the  a/ p ratios  ~10  provided  by  the 

low-altitude  measurements  discussed  above.  Several  parameters  of  the  model  were 

adjustable,  e.g.,  the  parameters  K and  t in  (1)  and  the  postulated  degree  of  mutual  ionic 

"thermalization"  (see  above)  at  L = 7,  and  it  was  not  very  difficult  to  achieve  a fit  of  the 

equatorial  theory  to  the  off-equatorial  data.  An  example  is  shown  in  Figure  12,  for 

92 

which  the  data  points  were  provided  by  Krimigis  . The  data  points  thus  correspond  to  L 
values  beteen  2.5  and  3.5  and  B values  in  the  range  0.18-0.22  G;  better  resolution  than 
this  was  unavailable  at  the  time. 

Except  for  the  fact  that  prohibitive  computing  costs  precluded  analysis  of  more 

35 

than  a single  value  of  J,  extension  of  the  above-described  work  to  particles  of  constant 
o 

(but  nonvanishing)  J /M  would  have  entailed  no  insurmountable  difficulty.  It  would  have 
been  necessary  to  have  averaged  the  atmospheric  density  over  bounce  motion  as  well  as 

drift  in  order  to  determine  the  Coulomb  and  charge-exchange  loss  rates,  but  this  had 

43  134 

already  been  done  . One  would  have  to  keep  track  of  the  variation  of  the  equatorial 
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FLUX,  ARBITRARY  UNITS 


Figure  11.  Relative  helium-ion  and  proton  flux  profiles  predicted  for  K/A  = 0.5 
MeV/nucleon  by  Cornwall'*'’  for  K = 10,  t = 20  min,  Ta  = 2T^  at  L = 7,  and  a doubled 


charge-exchange  loss  rate. 


compiled  by  Krimigis  from  off-equatorial  observations. 


2 

pitch  angle  qq  and  kinetic  energy  E with  L at  fixed  J /M,  and  one  would  need  to  take 

57  134 

account  of  the  variation  ’ of  with  °0  (see  Figure  13),  as  well  as  with  E and  L. 

(e) 

It  turns  out  (as  was  noted  above)  that  D ^ depends  on  E/A,  Z/A,  and  oQ  only  through  the 

drift  frequency  u>d/2tt,  whereas  depends  on  oQ  through  a factor  that  is  associated 

54  91  134 

with  drift-shell  tracing  ’ ’ . The  extension  of  Cornwall's  numerical  analysis  to 

o 

particles  with  nonvanishing  J /M  should  be  given  top  priority  in  the  theoretical  study  of 
ionic  radiation  belts. 

The  difference  in  equatorial  pitch-angle  anisotropy  between  protons  and  alpha 
particles  has  been  discussed  above  in  the  context  of  mutual  ionic  "thermalization"  at  the 
source.  However,  one  must  also  take  account  of  differences  in  Coulomb  and  charge- 
exchange  loss  rates  (recall  Figure  10)  and  of  the  considerations  outlined  in  the  paragraph 
above  in  order  to  obtain  quantitative  results.  Since  was  chosen  to  be  very  small  in 

Cornwall's  numerical  work'*'’,  and  since  is  rather  insensitive  to  a at  fixed  E,  it 
might  seem  to  make  little  difference  that  the  particles  were  assigned  a Q = 90°  there. 
The  difficulty  arises  from  the  systematic  variation  of  with  L,  which  impacts  not  only 
the  variation  of  E with  L but  also  (as  a consequence  of  this)  the  variation134  of  with 

o 

L at  fixed  J /M  (an  energy-independent  invariant,  which  therefore  is  unaffected  by 
Coulomb  loss).  It  may  well  be  true  that  is  negligible  (see  Figure  13)  for  particles 

that  mirror  far  off  the  equator,  but  his  remains  to  be  seen. 

The  point  of  the  rather  limited  comparisons  made  here  between  heavy-ion  theory 
and  the  observations  that  have  so  far  occurred  is  not  to  exult  in  the  fits  that  have  been 
achieved,  but  rather  to  emphasize  the  prospect  of  settling  decisively  and  quantitatively 
the  uncertainties  that  abound  concerning  magnetospheric  processes.  There  are  just  too 
many  possibilities  for  adjusting  parameters  and  too  many  effects  that  trade  off  against 
each  other  to  do  this  by  observation  of  protons  alone.  The  great  sensitivity  of  many  such 
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Figure  13.  Quantitative  representation  ’ of  dependence  of  radial  diffusion  co- 
efficients and  on  equatorial  pitch  angle  and  charge  state.  The  relative 

displacement  of  right  and  left  ordinates  remains  to  be  determined  empirically  for  given 
E/A  and  L. 
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effects  to  Z and  A allows  them  to  be  separated  clearly.  There  is  a continuing  need  for 
improved  observations  in  this  context,  e.g.,  to  obtain  data  on  helium  ions  having  E/A 
~ 10-15  keV/nucleon  at  L > 5.  While  we  have  emphasized  helium  ions  here,  the  same 

CO 

general  picture  applies  to  other  ions  as  well.  Fritz  and  Wilken  have  observed  heavy 

ions  at  ATS-6  in  synchronous  orbit,  ions  which  they  tentatively  identify  as  oxygen  of 

149 

solar-wind  origin.  Spjeldvik  and  Fritz  have  calculated  the  effects  of  magnetospheric 
radial  diffusion  and  charge  exchange  on  solar-wind  oxygen  ions  in  various  charge  states, 
in  analogy  with  the  work  of  Cornwall35  on  H+,  He+,  and  He++. 


IV.  MULTI-ION  PLASMA  PHYSICS 


The  presence  of  several  distinct  ionic  species  in  a collisionless  cold  plasma 

expresses  itself  by  introducing  new  cutoffs  (n  = 0 ) and  resonances  (n  = °°)  in  the 

refractive  index  n,  as  determined  by  magneto-ionic  theory'*®’*^4.  In  addition  to  the 

expected  ion-cyclotron  resonances  and  ion-ion  cutoffs  for  electromagnetic  waves 

propagating  parallel  to  the  magnetic  field  EL  one  obtains  various  ion-electron  and  ion-ion 

hybrid  resonances  for  oblique  propagation  . A major  historical  purpose  of  magneto-ionic 

124 

theory  has  been  to  study  the  propagation  of  radio  waves  through  the  ionosphere 

Although  the  proton  gyrofrequency  falls  well  below  the  VLF  (very-low-frequencv; 

u)  /2it  = 3-30  kHz)  band  at  ionospheric  altitudes,  one  may  expect  to  see  (ion-electron) 

hybrid-resonance  effects  there  at  VLF.  The  more  exotic  ion-ion  effects  might  occur  at 

frequencies  below  the  VLF  band,  e.g.,  at  3-3000  Hz  (commonly  known  as  the  ELF,  or 

extremely-low-frequency,  band)  or  at  frequencies  w/2tt  s 3 Hz  (commonly  known  as  the 

ULF,  or  ultra-low-frequency,  band).  Thus,  one  might  expect  ionospheric  H+,  He+,  and  0+ 

to  affect  the  downward  propagation  of  ULF  communication  signals  transmitted  from  the 

magnetosphere  and  intended  for  submarine  reception.  The  same  ions  might  affect  the 

82 

propagation  to  earth  of  natural  geomagnetic  pusations  generated  in  the  magneto- 
sphere. 

Magneto-ionic  resonances  (n  = « ) and  cutoffs  (n  = 0 ) are  often  associated  with 
band  limits  on  radio  noise  observed  in  space.  For  example,  the  VLF  receiver  aboard  the 
Alouette  1 spacecraft  observed  emissions  that  were  characterized  by  a sharply  defined 
lower  band  limit  identified  with  the  lower  hybrid  resonance1®’^4, Satellite 
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observations  of  chorus  and  ELF  similarly  show  a sharp  lower  band  limit,  but  this  has  been 

7 9 1 17 

identified  with  a two-ion  cutoff  in  the  extraordinary  wave  mode  ’ . Observations  of 

VLF  electric  and  magnetic  fields  with  the  lnjun-5  satellite  revealed  a host  of  sharply 
defined  noise  bands,  both  broad  and  narrow,  that  were  clearly  related  to  several  of  the 
various  cutoffs  and  resonances  obtained  in  multicomponent  magneto-ionic  theory7"*. 
Such  observations  yield  inferences  concerning  the  relative  abundances  of  various  ions  in 
the  plasma,  since  the  resonance  and  cutoff  frequencies  themselves  depend  upon  the  ionic 
masses  and  fractional  concentrations11,140. 

Notwithstanding  the  successes  noted  above,  magneto-ionic  theory  constitutes  only 
one  aspect  of  magnetospheric  heavy-ion  plasma  research.  The  main  thrust  in  recent 
years  (and  for  the  foreseeable  future)  has  been  toward  understanding  the  effects  of  heavy 
ions  on  magnetospheric  plasma  instabilities.  Such  effects  may  or  may  not  entail  a 
resonance  between  a wave  (w,J<)  and  a heavy  ion.  A wave-particle  resonance  would  be 
characterized  by  attainment  of  the  condition 


k||V(|  = ; ( = 0,  ±1,  ±2,  . . . 


(2) 


in  the  case  of  a nonrelativistic  particle  of  velocity  j/^and  (signed)  gyrofrequency  S2 j/2-rr  . 
The  subscript  j identifies  the  species  and  charge  state  of  the  ion;  the  subscript  n 
denotes  the  projection  of  the  associated  vector  in  the  direction  of  the  magnetic  field  B. 
The  case  f = 0 in  (2)  is  known  as  the  Landau  resonance  because  of  its  historic  role  in  the 


theory  of  damped  electrostatic  oscillations  in  an  electron  plasma 
correspond  to  the  various  orders  of  cyclotron  resonance. 


97 


The  cases  ( t 0 
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For  the  case  of  an  electromagnetic  wave  propagating  (i.e.,  having  its  ^0  in  the 
direction  of  B,  only  the  principal  cyclotron  resonance  (f  = ±1,  depending  on  the 

M* A 

direction  of  circular  polarization  of  the  wave)  is  important.  For  an  electrostatic  wave 
with  k parallel  to  J3,  only  the  Landau  ( t = 0)  resonance  is  important.  The  other 
resonances  in  (2)  occur  only  for  oblique  or  perpendicular  propagation  (kxB^O).  It 
follows  from  (2)  that  there  is  a certain  minimum  energy 


E*  = m.  v 2/2 
t ) ti 


(m./2k2  )(M.  - w)2 


(3) 


required  for  a particle  of  mass  m^  to  be  resonant  (in  order  t ) with  a wave  of  frequency 
uj/2tt  . In  the  case  of  a hot  plasma  there  will  normally  be  some  such  resonant  particles 
for  every  wave  (except  that  admissible  values  of  t are  severely  restricted  at  parallel 
propagation,  as  is  noted  above).  In  the  limit  of  a cold  plasma,  resonance  can  occur  only 
for  w = ± or  for  k(|  = ®,  i.e.,  under  conditions  corresponding  to  cold-plasma 
resonances  (n  = “> ) in  magneto-ionic  theory. 

Resonance  in  the  sense  of  (2)  leads  to  an  exchange  of  energy  between  the  wave 

and  the  plasma.  In  a hot  plasma  this  can  lead  either  to  growth  (instability)  or  to  damping 

of  the  wave.  Charged-particle  beams  and  other  velocity-space  anisotropies  are  major 

sources  of  free  energy  for  producing  instabilities,  as  are  gradients  in  the  spatial  distribu- 
77  78  110 

tion  of  plasma  Plasmas  satisfying  the  Maxwell-Boltzmann  distribution  are 

inherently  stable  and  cannot  give  rise  to  wave  growth,  i.e.,  to  an  amplification  of  the 
energy  content  of  a wave.  Gurnett  has  identified  electromagnetic  noise  bands  (in  the 
Hawkeye-1  data)  that  seem  to  have  arisen  through  resonant  amplification  by  anisotropic 
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populations  of  energetic  H+,  He++,  and  heavier  ions  in  the  equatorial  magnetosphere. 
86 

Kindel  and  Kennel  have  analyzed  the  conditions  for  instability  of  the  electrostatic  ion- 

acoustic  and  ion-cyclotron  wave  modes  in  a realistic  model  of  the  (multicomponent) 

topside  ionosphere.  The  free  energy  in  their  analysis  was  provided  by  a current  of 

precipitating  electrons,  but  resonances  between  the  waves  and  the  heavy  ions  were 

essential  to  the  results  obtained.  In  fact,  such  ion-cyclotron  waves  have  recently  been 

invoked  to  account  for  the  inferred  preferential  heating  of  upward-flowing  auroral 
i fin  i fi i 

ions  beams  ’ in  the  temperature  component  perpendicular  (rather  than  parallel)  to 
B. 

Heavy  ions  need  not  be  resonant  with  a wave  in  order  to  affect  the  growth  rate. 

o 

The  presence  of  cold  heavy  ions  (j  = h),  for  example,  can  affect  the  value  of  belonging 
to  a wave  that  is  resonant  only  with  hot  protons,  i.e.,  with  j =p  in  (3).  The  case  of  an 
electromagnetic  ion-cyclotron  wave  with  k parallel  to  B yields  a spatial  growth  rate 

- 1m  k|( « (2"2Npqp2/n2u>2)(2irmp*Tp  f* 

x[(TAP/Tp)(Op-w)  -np]sgnk|( 

x exp[-(mpc2/2*Tpn2u>2)(Slp  - w)2]  (4) 

42 

under  these  circumstances  , where  mp  , qp,  and  Np  are  the  proton  mass,  charge,  and 
number  density,  respectively.  The  proton  temperature  is  assumed  anisotropic  with 
respect  to  B^  such  that  /rTp  > tfTp  , where  k is  Boltzmann's  constant.  The  presence  of 
a cold  heavy-ionic  constituent  (j  = h)  with  its  complement  of  electrons  alters  Im  kj(  by 
changing  n2: 


W 

u 
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A (n2)  a (4 it  Nhqhc/B)  («  - w)'1  . 


(5) 


Thus,  if  the  argument  of  the  exponential  function  in  (4)  exceeds  unity  in  absolute 

iO  40  101 

value*  ’ ’ , then  the  addition  of  cold  heavy-ionic  plasma  will  enhance  the  spatial 

36  42  102 

growth  rate  ’ ’ for  waves  having  w < ft  ^ but  diminish  the  spatial  growth  rate  for 

u>nh- 

It  is  more  difficult  than  this  to  estimate  the  effects  of  "cold"  plasma  on  electro- 
static wave  modes,  partly  because  oblique  propagation  is  required  in  order  to  generate  a 
situation  of  physical  interest.  The  search  for  maximum  spatial  growth  in  this  case 

reduces  in  part  to  a numerical  search  for  minima  in  the  group  velocity  9u;/ak  and 

3 4 

ultimately  for  nonconvective  instability  ’ . Since  the  dispersion  relation  is  contingent  on 
plasma  temperature  here  to  a greater  extent  than  in  the  electromagnetic  case,  it  is 
important  not  to  set  T^  = 0 for  any  of  the  constituents.  However,  a mixture  of  hot  and 
"cold"  plasmas  seems  able  to  produce  nonconvective  instability  where  the  hot 
(anisotropic)  plasma  alone  would  produce  only  convective  instability.  The  foregoing 
comments  are  abstracted  from  a study  of  the  effects  of  hot  and  "cold"  electrons  on 
upper-hybrid  waves  ’ . Similar  work  has  been  done  recently  on  the  effects  of  hot  and 
"cold"  ions  of  a single  species  on  the  electrostatic  ion-cyclotron  loss-cone 
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instability  . This  instability  should  operate  outside  the  plasmasphere,  and  especially 
in  the  auroral  region.  Since  it  entails  oblique  propagation  (nonvanishing  k (|  ),  it  is  subject 
to  Landau  damping  by  "cold"  electrons  having  thermal  velocities  comparable  to  those  of 
the  hot  ions5.  Thus,  the  electrostatic  ion-cyclotron  instability  probably  does  not  operate 
inside  the  plasmasphere,  where  the  corresponding  electromagnetic  instability  tends  to  be 


enhanced  by  virtue  of  (5).  In  view  of  the  importance  of  electrostatic  cyclotron 
instabilities  outside  the  plasmasphere,  however,  it  would  be  of  great  interest  to  extend 
past  analyses  to  the  case  in  which  hot  and  "cold"  ions  of  several  species  are  present 

simultaneously.  Certain  aspects  of  the  Post-Rosenbluth  electrostatic  loss-cone  instabil- 
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ity  ’ have  been  analyzed  in  this  context  , but  much  more  work  remains  to  be  done 

on  this  important  subject.  A brief  review  of  density-sensitive  instabilities  in  general  is 
38 

given  by  Cornwall  . 

Loss-cone  anisotropy  is  a natural  consequence  of  mirror  geometry.  However,  one 

does  not  usually  find  an  otherwise  isotropic  pitch-angle  distribution  truncated  abruptly  at 

the  critical  (loss-cone)  angle  corresponding  to  particle  access  to  the  dense  atmosphere. 

Observed  distributions  more  nearly  resemble  the  bi-Maxwellian  form  used  in  (4),  insofar 

as  angular  distribution  is  concerned.  This  configuration  would  be  a natural  consequence 

of  pitch-angle  diffusion  in  the  presence  of  a loss  cone,  since  the  anisotropy  would  then  be 

distributed  throughout  the  distribution,  which  can  be  decomposed  into  eigenfunctions 

1 26 

that  vanish  at  the  edge  of  the  loss  cone  . 

Additional  sources  of  anisotropy  exist  for  both  light  and  heavy  ions.  One  of  these 

is  radial  diffusion  at  fixed  M and  J (the  first  two  adiabatic  invariants).  The  fractional 

energy  gain  associated  with  inward  transport  in  L from  an  external  source  varies 

inversely  with  mirror  latitude.  Thus,  if  one  examines  the  equatorial  pitch-angle 

distribution  at  fixed  E and  L,  the  particles  with  the  smaller  pitch  angles  must  have 

originated  in  the  higher-energy  (and  presumably  more  impoverished)  portions  of  the 
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source  spectrum  ’ . The  other  source  of  ion  anisotropy  is  charge  exchange. 

Charged  particles  experience  a bounce-averaged  atmospheric  density  that  increases  with 
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mippor  latitude  , and  so  the  particles  with  the  smaller  equatorial  pitch  angles  are  the 
most  rapidly  depleted  by  charge  exchange. 

Such  anisotropies  are  subject  to  limitation  by  the  plasma  instabilities  that  they 

produce  (see  above).  An  electromagnetic  proton-cyclotron  instability  resulting  from 

39 

proton  charge  exchange  has  been  analyzed  by  Cornwall  . The  consequent  pitch-angle 

diffusion  leads  to  a more  rapid  loss  of  protons,  and  to  a less  rapid  growth  of  proton 

anisotropy  with  time,  than  charge  exchange  alone  would  produce  during  the  recovery 

phase  of  a magnetic  storm.  Indeed,  such  ion-cyclotron  waves  have  been  observed  at 

large  amplitude  (0.4- 6.0  y at  1-30  Hz)  in  the  equatorial  magnetosphere  (L  ~ 3-5)  in 

153 

association  with  enhanced  anisotropic  ion  fluxes  on  several  occasions  . However,  the 

combination  of  proton  charge  exchange  and  pitch-angle  diffusion  ought  to  produce  a 

more  rapid  loss  of  particles  than  is  inferred  from  the  observed  decay  of  the  ring  current 

itself,  unless  the  recovery-phase  ring  current  consists  mainly  of  heavier  ions  such  as 

helium101**'*®.  It  even  seems  that  charge  exchange  alone  is  more  than  sufficient  to 

account  for  the  observed  decay  of  a ring  current  assumed  to  consist  entirely  of  protons 
146 

and  electrons  . Helium  ions  have  a smaller  charge-exchange  cross  section  for  neutra- 
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lization  by  the  atomic-hydrogen  atmosphere  at  E < 50  keV  than  do  protons  ’’  ; 

thus,  helium  may  have  become  the  dominant  ring-current  constituent  by  default  during 
the  recovery  phase  (see  Figure  10,  in  Secton  III). 

The  prospect  of  a helium-dominated  ring  current  during  recovery  phase  at 
E < 50  keV  certainly  requires  a new  analysis  of  the  ion-cyclotron  instabilities  that  might 
occur.  During  the  transition  from  hydrogen  dominance  to  helium  dominance,  for 
example,  one  must  take  account  of  both  ionic  species  in  the  hot-plasma  dispersion 
relation.  Moreover,  the  resulting  ion-cyclotron  waves  can  serve  as  a medium  for  energy 
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exchange  between  the  ionic  species.  The  cited  argument*®*  for  a helium -dominated 
recovery  phase  ring  current  ignores  the  role  of  radial  diffusion,  which  presumably 
persists  even  while  the  ring  current  is  decaying.  However,  radial  diffusion  seems  to 
favor  proton  dominance  only  for  drift  periods  «.  2 hr,  i.e.,  for  E » 30(Z/L)  keV  (see 
Figure  12  in  Section  III).  It  follows  from  (1)  that  is  substantially  independent  of 

species  for  ring-current  energies  (E  < 50  keV)  that  are  of  interest  in  the  present 
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context  . Observations  of  precipitating  ions  during  recovery  phase  have  revealed 

little  helium  but  much  hydrogen  and  oxygen.  Direct  observations  of  the  equatorial  ring 

current  with  ion-identifying  spectrometers  are  not  yet  available.  However,  the  question 

of  helium  dominance  during  recovery  phase  might  well  be  resolved"**  by  examining  the 

spectrum  of  equatorial  ions  observed80, ***** 14,129  at  very  iow  altitudes,  earthward  of 

the  inner  zone.  These  ions  form  a partial  radiation  belt  that  is  fed  by  the  stripping  of 

energetic  neutral  atoms  that  were  created  at  high  altitudes  ( » 20000  km)  by  the  charge 

exchange  of  ring-current  ions.  Identification  of  species  is  made  possible  by  the  fact  that 

the  product  of  cross  sections  for  the  two  steps  in  the  process  depends  much  differently 

on  energy  for  the  two  major  ring-current  ions*1*. 

When  considering  the  effects  of  heavy  ions  on  plasma  instabilities,  one  should  not 

on 

discount  the  effects  of  the  same  instabilities  on  the  heavy  ions  themselves.  Cladis  , for 
example,  points  out  that  electromagnetic  proton-cyclotron  waves  destabilized  by 
magnetospheric  proton  anisotropy  may  energize  ambient  oxygen  ions  at  low  altitudes. 
He  invokes  the  phase-trapping  of  particles  by  such  waves  in  the  inhomogeneous  medium 
to  account  for  energization  beyond  that  possible  in  quasilinear  theory.  This  is 
reminiscent  of  an  analogous  mechanism  proposed  by  Swift*80  for  electron  energization 
by  downward-propagating  electrostatic  waves,  a mechanism  that  has  since  been 
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investigated  numerically  and  analytically  . Brice  and  Lucas  have  described  the 

heating  of  natural  helium  and  oxygen  ions  (and  of  artifically  injected  lithium  ions)  in  the 

ring  current  as  a natural  consequence  of  instability  enhancements.  Moreover,  the 

enhancement  of  such  instabilities  may  lead  to  the  nonlinear  coupling  of  large-amplitude 

waves,  (Wj.Jtj)  and  (wj.kj),  to  produce  disturbances  at  ( Wj  ± u,2’ Jil  ± that  may 

interact  resonantly  with  portions  of  the  particle  distribution  that  would  have  been 

20  22 

accessible  to  neither  wave  in  the  linear  theory  ' . Various  quasilinear  and  nonlinear 

41 

aspects  of  the  mechanism  whereby  unstable  electromagnetic  proton-cyclotron  waves 
heat  ambient  electrons,  which  conduct  their  added  energy  down  the  field  line  and  thereby 
excite  atomic  oxygen  to  produce  a SAR-arc  (see  Section  I),  have  been  discussed  by 
Galeev65. 

In  concluding  this  long  section  on  heavy-ion  plasma  physics,  we  digress  to  make  a 
more  general  remark  about  space  plasma  physics.  To  an  outside  observer  it  often  seems 
that  plasma  theorists  concoct  a plasma  instability  to  account  for  every  observation,  and 
that  we  search  for  an  observation  to  support  every  instability  that  theory  predicts. 
There  is  some  truth  in  this  accusation.  Perhaps  the  role  of  plasma  instabilities  as  such  is 
overstated.  Let  us  consider  what  the  consequence  of  a plasma  instability  might  be:  a 
spectrum  of  turbulent  fluctuations  that  can  interact  with  the  distribution  of  charged 
particles.  On  the  other  hand,  what  can  be  said  of  a stable  plasma  containing  substantial 
free  energy?  One  can  say  that  it  is  characterized  by  a spectrum  of  fluctuations  that 
interact  with  the  distribution  of  charged  particles,  a spectrum  that  becomes  more  and 
more  intense  as  the  parameters  of  the  system  are  adjusted  toward  a state  of  marginal 
instability In  terms  of  the  physical  consequences,  there  is  a continuous  transition 
between  stability  and  instability.  Only  the  mathematical  description  is  dichotomous. 
Let  us  proceed  with  ths  perspective  in  mind. 
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V.  ACTIVE  PLASMA-INJECTION  EXPERIMENTS 


As  has  been  discussed  in  the  previous  section,  a magnetospheric  plasma  that 
contains  more  than  one  ionic  species  in  substantial  abundance  may  have  quite  different 
properties  from  a plasma  without  the  heavy  ions.  These  differences  extend  to  both 
wave-propagation  and  instability  characteristics.  Some  of  the  properties,  usually  those 
concerning  emission  at  the  very  lowest  frequencies  (w  « ),  depend  strongly  on  the 

species  of  the  heavy  ion  that  is  present.  Other  properties,  notably  those  that  affect 
analogous  electron  emissions  (see  below),  depend  mostly  on  the  total  charge  density  of 
the  added  ions.  All  active  experiments  presently  contemplated  depend  on  photo- 
ionization to  make  plasma  from  injected  neutral  particles;  thus,  the  only  substances  that 
are  useful  in  this  context  are  those  having  photo-ionization  times  < 1 hr  in  sunlight. 
Hydrogen  is  excluded  by  this  consideration,  but  barium,  cesium,  lithium,  and  perhaps  a 
few  other  substances  are  practical  for  artificial  plasma  injection. 

The  use  of  artificial  plasma  injection  as  a technique  for  enhancing  magnetospheric 

18  19 

plasma  instabilities  was  first  proposed  by  Brice  ’ in  the  context  of  whistler-mode 

85 

(electron-cyclotron)  waves.  The  whistler-mode  instability  seems  to  be  driven  by  an 
electron-'*  yjrature  anisotropy  (Tx  > T(|  ),  and  the  spatial  growth  rate  (-Imk()  is 
that  given  by  analogy  with  (4): 

- Im  k„  « (2  TT2Neq^/n2u>2)  (2v  me#(T*  f * 

X [(TZ/T^MlnJ-w)  - | nel]sgnk(| 
x exp[-(mec2/2xT||en2w2)(|ne|- w)2],  (6) 
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where  the  subscript  e refers  to  the  hot  electrons.  The  presence  of  a cold  heavy-ionic 

O 1 

plasma  constituent  (j  = h)  with  its  complement  of  cold  electron  would  alter  Im  k|(  in 

o 

(6)  by  changing  n : 


A(n4)  * (4ir  Nhqh/me)(|  nj-  u>)_1(u>  * n/1  |qe|. 


%-l 


(7) 


2 

Thus,  the  enhancement  of  - Im  k|:  via  A(n  ) is  independent  of  ionic  species  for  ft. 


[0  6 *1 

1 - (T|(  /TL)J|fte)  and  contingent  only  upon  the  charge  density  Nhqh,  provided 


that  the  argument  of  the  exponential  in  (6)  exceeds  unity  in  absolute  value49. 

The  amount  of  plasma  needed  to  conduct  an  active  experiment  is  not  large. 

Typical  barium  releases  inject  a fraction  of  a kilogram  of  Ba,  which  is  photo-ionized  in 

'-20  sec.  At  an  expansion  velocity  ~ 1 km/sec,  an  ionospheric  plasma  cloud  would  thus 

attain  a nominal  radius  ~ 20  km  across  B,  but  it  might  ultimately  attain  a length 

> 103  km  along  JB.  Half  a kilogram  of  Ba+  would  fill  this  volume  to  a density 
—3 

~ 2000  cm  , and  even  higher  densities  would  have  been  attained  before  the  barium 

cloud  reached  its  full  length.  There  are  three  methods  of  Ba+  injection:  (1)  rocket-borne 

canisters  of  barium  exploded  in  the  ionosphere  to  produce  a cloud  that  is  initially 
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isotropic  and  confined  to  the  ionosphere  ’ ’ ; rocket-borne  shaped-charge  releases, 

which  create  a cloud  that  can  extend  along  an  entire  field  line*63  *99.  anq  (3)  high- 
altitude  releases  in  the  magnetosphere  . There  is  a great  difference  between  the 
physics  of  ionospheric  Ba+  releases  and  the  physics  of  high-altitude  releases.  Ionospheric 
release  yields  a low  value  (« 1)  for  p , which  is  defined  as  the  ratio  of  plasma  pressure  to 
B /8t r.  in  this  case  there  is  only  a minimal  distortion  of  the  earth's  local  magnetic  field 
B.  There  may  or  may  not  be  a major  distortion  of  the  local  electric  convection  field  £. 
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This  would  depend  on  the  resistance  of  the  electric  circuit  that  links  the  plasma  cloud 
with  the  ionosphere.  If  there  is  no  major  distortion  of  £,  then  the  barium  cloud  drifts^4 

O 

with  the  ambient  convection  velocity  (c/B  ) £ x B.  Moreover,  there  occur  drift-wave 

«*»  ****** 

instabilities  that  lead  to  an  easily  observed  striation  of  the  plasma  cloud1®®’1®4  in  the 
presence  of  collisions  between  the  Ba+  ions  and  the  neutral  atmosphere.  A barium 
release  at  high  altitude  (e.g.,  several  earth  radii  out  in  the  magnetosphere)  does  not 
entail  collisions  between  the  Ba+  ions  and  the  neutral  atmosphere,  but  it  does  entail  a 
high-beta  ( p » 1)  condition.  The  result  is  a severe  distortion  of  the  local  B field,  such 
that  both  B and  £ are  initially  excluded  from  the  expanding  plasma  cloud,  which  is 
highly  conducting.  Eventually  the  earth's  magnetic  field  diffuses  back  into  the  plasma 
cloud  and  the  electric  convection  field  penetrates  the  cloud  as  well.  Meanwhile,  the 
cloud  has  transferred  its  excess  momentum  (over  that  which  it  would  have  contained  if 
traveling  at  the  convection  velocity)  to  the  ambient  plasma.  For  plasma  clouds  such  as 
those  released  from  HEOS  1 at  an  altitude ~70000  km,  however,  this  transfer  takes  place 
on  a time  scale  of  hours  . Scholer  and  Pilipp  have  discussed  the  motion  to  be 
expected  from  such  magnetospheric  barium  releases. 

It  has  not  been  possible  so  far  to  draw  any  definite  conclusions  about  the  influence 
of  such  high-altitude  releases  on  the  magnetosphere  itself.  The  barium  clouds  have  such 
small  transverse  dimensions  that  it  is  difficult  to  locate  (by  remote  means)  the  tiny 
volume  of  magnetosphere  that  would  be  affected.  (However,  Wescott  et  al.  have 
discussed  a substorm  that  may  possibly  have  been  triggered  by  a shaped-charge  barium 
release.)  On  the  other  hand,  there  is  some  evidence  for  increased  wave  emission  and 
increased  energetic-electron  flux  following  ionospheric  barium  releases,  which  often  find 
the  instrumented  rocket  in  the  middle  of  the  plasma  cloud.  Such  effects  have  been 
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reported  by  Kelley  et  al.  and  by  K8hn  and  Page  ; the  effects  were  too  complex  to  be 

90 

described  in  detail  here.  More  recently,  Koons  and  Pongratz  have  reported  the  direct 
excitation  of  electrostatic  ion-cyclotron  waves,  in  the  ambient  ionospheric  plasma  as 
well  as  in  the  artificially  injected  plasma,  by  the  impulse  associated  with  a shaped- 
charge  Ba  release  at  450  km  altitude. 

Although  it  is  very  much  worthwhile  to  continue  experiments  that  involve  the 
release  of  barium,  one  should  recognize  that  barium  has  certain  limitations  for  use  in 
magnetospheric  research.  The  small  transverse  dimension  has  been  noted  above, 
together  with  the  high- (3  conditions  that  this  often  implies.  Moreover,  barium  plasma  is 
unsuitable  for  experiments  designed  to  enhance  electromagnetic  ion-cyclotron  instabil- 
ities by  means  of  (5).  To  create  a meaningful  enhancement  of  - Im  k , one  would  wish  to 
inject  cold  plasma  outside  the  plasmasphere  (where  there  is  little,  if  any,  cold  plasma  of 
natural  origin)  but  near  the  equator  (where  the  hot  plasma  is  concentrated  by  virtue  of 
the  inequality  TP /Tp  > 1).  The  Ba+  gyrofrequency  can  hardly  exceed  0.05  Hz  there, 
and  only  waves  below  this  frequency  can  have  their  growth  rates  enhanced  by  the 
addition  of  cold  plasma  . Except  perhaps  within  the  barium  cloud  itself,  the  wavelength 
corresponding  to  such  a low  frequency  will  likely  exceed  the  diameter  of  the 
magnetosphere.  It  is  difficult  to  imagine  observing  such  a "wave"  at  all,  not  to  mention 
its  "enhanced"  instability.  Moreoever,  the  growth  rate  at  such  low  frequencies  is  very 

small  compared  to  the  maximum  that  occurs  at  w = [1  - (TP  /TjP  )]  t in  the  absence 
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of  cold  plasma  ’ , where  0.8  < « < 1 if  P < 1.  Thus,  the  addition  of  cold  barium 

plasma  to  the  magnetosphere  may  well  enhance  electron-cyclotron  instabilities  via  (7), 
but  would  tend  to  suppress  the  ion-cyclotron  instabilities  that  are  (in  some  respects) 
more  interesting.  The  same  objection  holds  in  the  case  of  cesium. 


As  for  lighter  ions,  the  only  parent  neutral  with  a reasonably  short  photo- 
ionization lifetime  (~lhr)  is  lithium,  which  has  repeatedly  been  released  into  the 
ionosphere  over  the  years,  both  from  rockets  and  through  thermonuclear  explosions. 
Neutral  Li  is  readily  observed  as  a resonant  scatterer  of  sunlight,  but  the  Li  ion  is 

completely  invisible.  Lithium  clouds  injected  from  rockets  launched  at  twilight  have 
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been  observed  with  brightnesses  in  the  kilorayleigh  range;  1 kR  corresponds  to  10  photon 

emissions  per  second  by  the  atoms  contained  in  a column  whose  cross  section  (measured 

o 

in  the  plane  transverse  to  the  line  of  sight)  is  1 cm  . 

It  has  been  proposed  (but  not  yet  attempted)  to  make  large-scale  lithium  releases 
near  the  equator  outside  the  plasmasphere  as  a means  of  triggering  electromagnetic 
instabilities  driven  by  ring-current  ions.  The  naturally  occurring  plasma  there  is  of  low 
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density  (often  < 1 cm  ),  and  it  is  technologically  feasible  to  add  lithium  plasma  at  a 
density  > 1 cm  over  a large  region  of  space  (see  below).  Readily  observed  effects 
should  occur  at  Li  densities  much  more  modest  than  1 cm  , but  the  enhancement  of 
-Im  k„  occurs  only  below  the  lithium  gyrofrequency  (w  « 0.14 12  p ) rather  than  at  the 
growth-rate  peak  in  the  absence  of  lithium  (w  ~ 0.4  12  p for  T^p  ~ 2T(|P  and  (3  ~ 1). 
For  example.  Figure  14  shows  the  effect  of  a 10%  increase  in  the  plasma  number 
density^2on  the  temporal  growth  rate  y = -(du>/dkj()  Im  kH  when  Li+  is  added  to  an 
already  unstable  hot  hydrogenic  plasma  having  = 2TP  and 

pp  = 8 ir  Np  kT^/B2  = 1.  The  result  is  a narrow  band  of  emissions  just  below  the  Li+ 
gyrofrequency,  where  there  had  been  none  before.  However,  the  growth  rate  for  this 
band  remains  much  smaller  than  the  peak  y in  the  absence  of  Li+,  and  one  would  need 
considerably  more  than  10%  lithium  plasma  (by  number  density)  to  make  the  peak  at 
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w a 0.14  0 dominant  over  that  at  w ~ 0.4  fi.  . These  considerations  have  led 
P P 
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Mflrk  to  advise  against  adding  lithium  to  an  already  unstable  hydrogenic  plasma 

outside  the  plasmasphere.  Of  course,  if  the  hot  ions  were  helium  (He+),  the  naturally 

occurring  emissions  would  peak  at  w ~ 0.1  Op  and  the  addition  of  even  small  amounts 

of  lithium  would  enhance  this  peak.  However,  the  arguments  cited  in  Section  IV  for  a 

helium-dominated  ring  current  are  contingent  on  charge  exchange,  which  decreases  in 

importance  with  increasing  L,  and  are  addressed  specifically  to  the  recovery  phase  of  a 

magnetic  storm.  One  cannot  rely  on  the  ring  current  to  be  helium-dominated  at  all  times 

and  in  all  places.  Far  more  dramatic  would  be  the  effect  of  adding  lithium  to  a stable 

hydrogenic  plasma  (e.g.,  having  < 1.2^  with  P|f~  1).  In  this  case  the  peak  value 

of  - Im  k(|  (which  probably  would  occur  at  w s 0.14  flp  ) is  likely  to  be  insufficient  (in 

the  absence  of  cold  plasma)  to  overcome  the  losses  that  occur  at  the  ends  of  the 

amplification  path.  The  addition  of  even  a modest  density  of  cold  Li+  should  enhance 

this  peak  value  of  - Im  k(|  enough  to  overcome  such  losses;  a plasma  having  roughly  equal 

number  densities  of  cold  Li+  and  hot  H+  has  been  found  to  be  unstable  for  anisotropies 
42 

that  were  quite  small  . 

Because  neutral  lithium  has  a photoionization  time  ~ 1 hr  and  a larger  thermal 

4 

velocity  than  similarly  injected  barium,  it  should  attain  a radius  ~ 10  km  transverse  to 

B.  A spherical  cloud  of  this  radius  would  require  ~ 150  kg  of  Li+  to  reach  a uniform 
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density  ~ 3 cm  . Such  a cloud,  however,  would  envelop  - 10  erg  (roughly  equivalent 

to  2.4  kilotons  of  TNT)  in  the  form  of  geomagnetically  trapped  (particle)  radiation,  some 

fraction  of  which  might  well  be  released  in  the  form  of  waves,  precipitating  particles, 

and  various  substorm-like  manifestations  of  the  expected  plasma  instabilities  if  the 

experiment  were  conducted. 
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In  fact,  injections  of  ions  lighter  than  barium  have  undoubtedly  occurred  already 

near  syncnronous  altitude  (~  35000  km),  when  rocket  engines  were  burned  to  maneuver  a 

119 

satellite  into  the  desired  circular  orbit.  Paulikas  has  looked  for  effects  associated 

with  such  injections  and  found  none,  but  probably  because  the  plasma  clouds  were  too 

small.  Kivelson  and  Russell  have  looked,  with  considerable  success,  for  such  effects 

associated  with  the  ’’injections"  of  cold  hydrogen  plasma  that  occur  naturally  during 

sporadic  erosions  of  the  plasmasphere  that  are  associated  with  changes  in  the  electric 
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convection  field.  Moreoever,  Kintner  and  Gnrnett  have  identified  from  their 
Hawkeye-1  data  several  episodes  of  electromagnetic  proton-cyclotron  wave  occurrence 
near  the  plasmapause,  often  just  inside  the  plasmasphere. 

In  some  respects  it  is  unfortunate  that  the  controlled  injection  of  H plasma  into 

the  magnetosphere  by  conventional  means  is  an  impossible  experiment.  Much  of  the 

motivation  for  active  plasma-injection  experiments  is  to  simulate  the  effects  of 

naturally  occuring  regions  of  high-density  magnetospheric  plasma.  It  was  originally 

33 

suggested  by  Cocke  and  Cornwall  that  the  electromagnetic  proton-cyclotron  mode 

could  be  made  unstable  by  a high  density  of  cold  plasma  in  the  magnetosphere.  This 

40 

picture  later  developed  into  a theory  of  storm-time  dynamics  and  of  SAR-arc 
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generation  . Some  results  ’ from  S (Explorer  45)  strongly  suggest  that  the 

equatorial  proton  ring  current  (E  ~ 1-200  keV)  interacts  with  the  plasmasphere  in 

essentially  the  manner  predicted  by  the  theory,  although  detailed  information  on 

154 

electromagnetic  cyclotron  waves  at  the  equator  is  sparse  . For  example,  Williams  and 
172 

Lyons  have  reconstructed  a hypothetical  plasma-density  profile  at  the  plasmapause 

Ifrom  their  particle  data  by  invoking  (3)  at  marginal  stability  for  1=1.  It  follows  from 
(4)  that  Im  k^  =0  for 


« = [i-(Tp  /tp  )]np  . 


(8) 


u>  = 
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The  refractive  index  n at  w = w can  be  calculated  as  if  in  a cold  plasma  (p  « 1). 

# 

Thus,  one  obtains  (for  w = w ) the  fact  that 

n2  s (ck|(  / w )2  = 1 + [l  -(«/Op  )]_1  (c/cA)2  , (9) 

n 

where  c^  is  the  Alfven  speed,  and  from  (3)  the  result  (assuming  n » 1)  that 
E*  « (m  p c 2/2)[(TP  /T  P ) - l]  "2  (TP  /T  ^ ) 

a (B2/8tt  Np)  [(TjP  /T  P ) - l]"2  (T  P /TP  ) , (10) 

where  Np  is  the  total  density  of  hot  and  cold  protons  (or  their  mass-equivalent)  in  the 

equatorial  plasma.  By  identifying  the  critical  energy  that  corresponds  to  a marked 
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change  in  the  pitch-angle  distribution  at  a given  L value,  Williams  and  Lyons  obtained 

~ 0 ♦ 

the  equatorial  density  profile  Np  =(Bq/8t:E1)  shown  in  Figure  15,  where  we  have 
eliminated  the  arbitrary  dependence  on  T p /T(|p  in  (10)  by  assuming  T^p  « 1.755  T p 
(this  being  reasonable,  if  not  self-consistent).  The  vertical  bar  represents  an  independent 
measure  of  Np  from  saturation  of  the  on-board  DC  electric-field  probe.  This  agrees 
well  with  the  profile  = (Bo/8ttEj)  , which,  moreover,  has  the  expected  shape. 

18  19 

The  early  ideas  of  Brice  ’ have  thus  undergone  considerable  elaboration  and 

o<*  JO  JI7  iq  1 An 

refinement  ’ ’ ’ during  the  past  several  years,  but  his  proposal  to  seed  the 

magnetosphere  with  cold  plasma  to  simulate  naturally  occurring  substorm  effects 


-59- 


5 


2 *.  17? 

Figure  15.  Plasma  densities  N ~ (Bq  /8itEj)  inferred  by  Williams  and  Lyons  from 


energy  sec  aQ  characteristic  of  transition  from  anisotropic  to  nearly  isotropic 


equatorial  pitch-angle  distribution  at  Explorer  45  (S  ).  Vertical  bar  denotes  independent 
density  estimate  based  on  observed  saturation  of  on-board  DC  electric-field  probe. 
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remains  as  attractive  as  ever.  Brice  did  not  specifically  take  into  account  the  damping 
effect  of  cold  heavy  ions  on  proton-driven  waves,  and  he  used  the  concept  of  a maximal 
stably  trapped  flux  rather  uncritically  to  predict  the  triggering  of  large-scale 
precipitation  of  energetic  particles  upon  the  injection  of  cold  plasma  into  the 
magnetosphere.  The  difficulty  of  destabilizing  proton-driven  instabilities  is  partially 

AC  JA 

overcome  by  using  lithium  in  place  of  barium  or  cesium  for  the  experiment00’  . Use  of 

lithium  rather  than  hydrogen  makes  it  impossible  to  simulate  naturally  occurring 

phenomena  exactly,  but  lithium  has  its  own  advantages.  For  a given  number  density  of 

+ 9 

ions,  lithium  is  seven  times  as  effective  as  H in  lowering  the  parameter  c^  in  (10),  and 

* 

thus  in  lowering  Moreover,  as  has  been  noted  above,  lithium  is  an  especially 

appropriate  additive  when  the  proton  anisotropy  is  small,  e.g.,  when  T^p  < 1.25  Tp  . In 

this  case  the  peak  growth  rate  in  the  absence  of  cold  plasma  might  well  have  occurred  at 
or  below  the  lithium  gyrofrequency,  a condition  that  would  overcome  the  objection  raised 
by  MHrk102  against  the  injection  of  Li+  into  a significantly  more  anisotropic 
(Tp  > 2TP  ) proton  ring  current.  The  performance  of  active  experiments  in  space  has 
a bright  future,  and  lithium-plasma  injection  into  the  magnetosphere  should  play  an 
important  role  in  that  future.  This  subject  will  be  addressed  again  in  Section  VII. 
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VI.  HEAVY-ION  PHYSICS  AT  JUPITER 


Needless  to  say,  the  study  of  heavy-ion  physics  in  other  planetary  magnetospheres 

is  even  less  advanced  than  it  is  in  the  earth's.  In  a recent  1200-page  comprehensive 

review  of  Jupiter  , for  example,  there  are  no  articles  specifically  addressed  to  the 

Jovian  analogue  of  any  of  the  previous  sections  of  the  present  paper.  Some  results  on 

energetic  alpha  particles  from  the  Pioneer-10  data145,158  have  been  reported  outside 

that  review.  For  energies  of  several  MeV/nucleon  the  a /p  flux  ratio  was  found  to  range 
-3  -2 

from  10  to  10  along  the  roughly  equatorial  orbit,  but  nothing  is  known  about  the  a/p 
ratio  at  lower  energies. 

The  possible  physical  mechanisms  that  might  involve  Jovian  heavy  ions  are  more 

elaborate  than  the  mechanisms  that  operate  in  the  earth's  magnetosphere.  In  addition  to 

ionospheric  and  solar-wind  sources,  one  must  consider  ions  derived  from  Jupiter's  major 

satellites,  which  orbit  well  inside  the  magnetosphere.  Many  more  ionic  species  are 

involved  at  Jupiter  than  at  the  earth,  and  some  of  these  (e.g.,  the  ions  of  He,  C,  O,  Ne, 

Na,  Mg,  S,  and  Fe)  can  be  used  for  distinguishing  sources,  acceleration  mechanisms,  and 
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transport  processes.  Neutral  sodium  has  been  observed  ’ around  the  orbit  of 

96  109  + 

Jupiter's  satellite  lo,  as  well  as  ionized  sulfur  ’ . The  clouds  of  Na  and  S appear  to 

form  a partial  torus  that  extends  eastward  and  westward  from  Io  and  encircles  the 
53  107  108 

planet  ’ ’ , thus  filling  a significant  region  of  Jupiter's  magnetosphere.  (Ionized 

sodium  and  neutral  sulfur  would  not  be  observable  by  earth-based  spectroscopy.)  No 
other  heavy  ion  has  been  definitely  identified  in  the  Jovian  magnetosphere.  The  other 
ions  listed  above  are  expected  on  the  basis  of  models  of  the  composition  of  Jupiter  and 
its  Galilean  satellites.  Io  is  also  known  to  have  a cloud  of  neutral  hydrogen  which 
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forms  a partial  torus,  and  Ganymede  (which  is  known  to  consist  largely  of  H.^O)  might 
theoretically  have  a similar  corona.  Such  clouds  of  neutral  hydrogen  could  play  an  im- 
portant role  in  removing  heavy  ions  by  charge  exchange,  just  as  in  the  earth's 
magnetosphere.  Aside  from  the  satellites,  there  are  essentially  no  known  sources  for 
low-energy  neutrals  in  the  Jovian  magnetosphere. 


One  outstanding  feature  of  Jupiter  is  the  rapid  rotation  rate  that  it  transmits  to 
the  magnetosphere,  seemingly  at  least  to  the  orbit  of  Ganymede  (r  w 1 5r j).  Indeed,  at 
particle  energies  < 8 MeV/charge  the  contribution  to  the  drift  rate  from  gradient- 
curvature  effects  is  smaller  than  the  corotation  rate  for  L s 15.  Particles  having  Z/Z 
« 8 MeV/charge  all  have  essentially  the  same  drift  frequency  and  (therefore)  the  sam 
radial-diffusion  coefficient  D^.  regardless  of  Z or  A;  tests  discussed  in  Section  III  to 
identify  the  radial-diffusion  process  are  inapplicable  at  Jupiter  because  of  this 


degeneracy.  Actually,  the  radial-diffusion  processes  represented  by  (1)  are  believed 

to  be  less  important  at  Jupiter  than  that  caused  by  the  fluctuating  electric  polarization 

23  44 

fields  associated  with  variable  neutral  winds  in  Jupiter's  ionosphere  ’ . 

The  rapid  corotation  of  Jupiter's  magnetosphere  drags  magnetospheric  plasma  and 

zenomagnetically  trapped  radiation  past  the  Galilean  satellites  so  as  to  produce  several 

possible  effects:  (1)  impact  of  corotating  energetic  particles  on  a satellite  may  liberate 

atoms  or  ions  into  the  magnetosphere;  (2)  corotation-induced  electrostatic  fields  may 

accelerate  thermal  plasma  (along  J3)  up  to  several  hundred  keV  at  Io,  or  up  to  90  keV  at 

Ganymede,  enabling  these  particles  to  liberate  atoms  or  ions  from  the  surface;  and  (3) 
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magnetospheric  ions  may  be  lost  not  only  through  direct  impact  on  a satellite,  but 

also  by  charge  exchange  as  they  drift  through  the  satellite's  corona  of  neutral  atoms. 


-64- 


Admittedly  little  is  known  about  the  rates  at  which  ions  or  atoms  can  be  sputtered  from 

a satellite  surface  by  a given  incident  particle  flux,  and  the  magnitude  of  the  incident 

particle  flux  itself  is  highly  uncertain.  It  has  been  pointed  out  that  if  Io  is  a good 
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electrical  conductor,  as  Piddington  and  Drake  have  proposed,  then  lo's  effective  cross 
section  for  impact  by  magnetospheric  particles  may  be  very  different  from  lo's 
geometric  cross  section.  Its  swath  is  found  to  be  enhanced  (over  the  geometric  value)  for 
electrons  having  E > 30  MeV.  A similar  effect  should  obtain  for  a satellite,  such  a 
Ganymede  (see  above),  that  is  a good  dielectric  (in  this  case  consisting  mostly  of  water). 

Charge-exchange  processes  have  implications  for  remote  sensing  of  satellite 

coronas  as  well  as  for  magnetospheric  dynamics.  A simple  measure  of  the  effectiveness 

of  charge  exchange  in  removing  ions  is  the  product  N cr  v t , where  N is  the  average 

number  density  of  neutral  atoms  in  the  corona,  t is  the  amount  of  time  that  the 

energetic  ion  spends  in  this  neutral  cloud,  and  is  the  charge-exchange  cross  section  at 

relative  ion-neutral  velocity  ^(which  is  essentially  the  ion  velocity  associated  with  the 

gyro  and  bounce  motion).  For  50-keV  protons  incident  on  atomic  hydrogen  (see  Figure 
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10)  one  obtains  uv  a 3X10  cm  /sec;  for  He  at  the  same  energy  one  obtains  av 
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« 1.5X10  cm  sec.  A magnetospheric  ion  might  require  t ~ 2 — 100X10  sec  to 

diffuse  across  a neutral  cloud  of  minor  diameter  2rj  at  Ganymede  (L  = 15)  and  perhaps 

15—200  times  as  long  through  a similar  cloud  at  lo  (L  = 6).  The  uncertainty  in  these  time 

estimates  reflects  to  some  extent  the  uncertainty  in  the  mechanism  primarily  responsible 

for  'ower  estimates  here  correspond  to  the  Brice-McDonough  mechanism 

(ionospheric  neutral  winds)  with  a recommended  coefficient 

~ 2.3  X 10”1  L/1  sec”1.  Significant  charge  exchange  of  H+  would  have  occurred  in  the 
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interim  if  N > 1/tvt,  i.e.,  if  N > 1 cm  around  lo's  orbit  or  if  N » 15  cm  around 


28  12-2 
Ganymede's.  Carlson  and  Judge  report  a column  content  ~ 3 x 10  cm  for  neutral 

hydrogen,  based  on  Lyman-o  measurements  from  Pioneer  10.  This  would  correspond  to  a 
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number  density  N ~ 200  cm  if  the  cloud  had  a diameter  comparable  to  that  of  Jupiter 
itself.  The  product  <tv  t decreases  drastically  with  increasing  energy,  but  significant 
charge-exchange  effects  might  still  occur  for  energies  up  to  several  hundred  keV  at 
Ganymede  and  several  MeV  at  Io. 

For  lack  of  space,  we  do  not  analyze  other  possible  injection,  loss,  and 
acceleration  mechanisms  for  heavy  Jovian  ions;  these  would  be  analogous  to  mechanisms 
that  operate  in  the  earth's  magnetosphere.  It  might  be  difficult  to  distinguish  between 
solar-wind  and  ionospheric  sources  for  trapped  Jovian  radiation,  since  Jupiter's 
atmosphere  should  have  a composition  rather  like  that  of  the  photosphere.  An 
interesting  exception  to  this  generalization  is  the  CNO  group,  whose  atoms  are 
presumably  tied  up  in  chemical  compounds  (e.g.,  H20,  NH3,  CH4>  C2H2)  in  Jupiter's 
atmosphere  and  have  low  elemental  abundances  there.  Apart  from  the  case  of  acetylene 
(CgHg),  however,  it  might  be  difficult  to  distinguish  the  ionized  atom  from  the  ionized 
molecule  in  terms  of  mass  or  charge. 

An  additional  note  concerning  sulfur  and  sodium  seems  in  order.  These  seem  to 
occur  in  a torus  surrounding  the  orbit  of  Io.  Strangely  enough,  sodium  and  sulfur  do  not 
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seem  to  co-exist  in  the  torus  ’ * ’ . Either  one  or  the  other  seems  to  be  present 

52 

at  any  given  time,  but  not  both.  It  has  been  suggested  that  the  transition  from  a sulfur 
torus  to  a sodium  torus  might  be  accompanied  by  a magnetohydrodynamic  (MHD) 
transformation  in  Jupiter's  magnetosphere,  i.e.,  that  Jupiter's  magnetosphere  is  like  the 
sun's  (with  strong  radial  outflow  of  plasma)  in  the  sulfur  phase  and  like  the  earth's  (with 
only  a weak  polar  wind  at  best)  in  the  sodium  phase.  Changes  in  the  Alfven  speed  cA 
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relative  to  the  equatorial  corotation  velocity  (13L  km/sec)  might  account  for  such  an 
52 

MHD  transformation  . 


Finally,  we  should  mention  that  the  study  of  Jupiter's  magnetosphere  teaches 
• 1 30 

important  lessons  about  the  earth's.  It  has  been  noted  that  the  radial-diffusion  coef- 
ficient D , corresponding  to  electrostatic  and/or  magnetic  impulses  of  the  sort 
described  by  (1),  should  vary  at  least  as  the  fifth  power  of  the  earth's  dipole  moment  Me> 

where  <t>  is  the  third  adiabatic  invariant.  (One  speaks  here  of  D rather  than  D. 
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because  L is  not  an  adiabatic  invariant  unless  Mg  is  constant  in  time.)  One  expects  the 
23 

Brice-McDonough  D to  have  a weaker  variation  with  Me>  perhaps  becoming  the 
dominant  radial-diffusion  mechanism  in  the  earth's  magnetosphere  at  epochs  of  unusually 
small  Me  (e.g.,  about  6000  years  ago,  when  Me  amounted  to  <50%  of  its  present  value). 
Thus,  a diffusion  mechanism  that  was  initially  contemplated  in  response  to  a puzzle 
presented  by  Jupiter  may  turn  out  after  all  to  have  played  an  essential  role  in  the  ancient 
history  of  the  earth. 


There  are  three  vital  observational  aspects  of  magnetospheric  heavy-ion  physics 
that  remain  unexplored  to  this  day:  (1)  as  was  mentioned  in  Section  1,  there  is  an  instru- 
mentation gap  for  ions  between  ~ 20  keV/charge  and  ~ 200  keV/nucleon,  and  no  direct 
charge-composition  measurements  are  available  in  this  range;  (2)  satellites  instrumented 
for  low-energy  ion-composition  measurements  have  not  traversed  the  region  5 £ L S 10 
near  the  equator,  a region  of  special  significance  in  view  of  the  very  sharply  peaked 
pitch-angle  distributions  of  protons  and  (especially)  alpha  particles  that  are  observed 
near  the  equator  at  lower  L values;  and  (3),  there  is  very  little  information  available 
about  the  distribution  of  heavy  ions  in  the  magnetospheres  of  other  planets.  A number  of 
recent  proposals  to  fill  these  observational  gaps  will  be  discussed  below. 
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Several  of  the  barium-injection  experiments  ’ ’ * cited  in  Section  V were 

essentially  passive  in  that  they  were  designed  to  trace  magnetic  field  lines  or  to  monitor 
drift  at  the  ambient  electric  convection  velocity  (c/B  ) £ x B.  More  experiments  of  this 
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type  are  planned  for  the  future.  Several  other  past  experiments  have  had  the  effect  of 
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modifying  localized  regions  of  the  magnetosphere  ’ ’ ’ , and  at  least  one  may 

have  triggered  a large-scale  magnetospheric  response.  Much  less  work  has  been  done, 
however,  with  unequivocally  active  plasma  injection  experiments  designed  to  modify  a 
large  region  of  the  magnetosphere.  While  some  of  the  proposed  experiments  are 
undoubtedly  too  expensive,  it  nevertheless  seems  quite  worthwhile  to  develop  at  least  a 
modest  program  for  large-scale  active  plasma  releases  in  the  magnetosphere. 

There  is,  of  course,  no  abrupt  transition  from  the  heavy-ion  physics  of  the  past  to 
that  of  the  future.  At  this  writing  the  Voyager  spacecraft  approaches  Jupiter  with 
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E*dE/dx  telescopes,  the  first  International  Sun-Earth  Explorer  (ISEE)  has  just  been 
launched  with  an  ion-composition  instrument  (E-dE/dx  telescope  plus  TOFD:  see  Section 
1),  and  the  GEOS  satellite  (with  ESA's)  had  preceded  ISEE.  These  satellites  undoubtedly 
will  provide  much  new  and  useful  data.  Unfortunately,  the  orbits  of  the  earth  satellites 
are  not  ideal:  although  they  reach  high  L values,  the  satellites  do  not  remain  near  the 
equator  for  any  significant  fraction  of  the  time.  Nor  do  they  fill  in  the  instrument  gap 
mentioned  above  and  in  Section  I. 

A number  of  groups  have  proposed  to  fill  most  of  the  instrument  gap;  a typical 
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solution  is  that  proposed  by  Gloeckler  : Incoming  ions  are  collimated,  and  those  of  a 
given  energy/charge  are  selected  by  an  electrostatic-deflection  analyzer.  These  ions 
pass  through  a TOFD  and  are  stopped  in  a detector  that  measures  their  total  energy. 
Energy,  mass,  and  charge  can  thus  be  measured  from  very  low  energies  up  to  about 
400  keV/charge.  The  instrument  is  capable  of  isotopic  separation  for  elements  as  heavy 
as  carbon  and  elemental  separation  up  to  iron.  The  upper  energy  limit  is  set  by  the  size 
of  the  TOFD  and  electrostatic-deflection  system.  An  instrument  of  this  type  will 
virtually  double  our  knowledge  about  low-energy  helium  ions  in  the  magnetosphere  by 
separating  He+  from  He++.  We  recall  from  Section  III  that  such  information  is  vital  in 
assessing  the  balance  between  radial-diffusion  and  charge-exchange  processes.  An 
unavoidable  drawback  of  the  above  instrument  is  its  rather  modest  upper  energy  limit. 
No  one  has  suggested  a way  of  measuring  charge  states  at  E/A  > 1 MeV/nucleon  with 
instruments  of  reasonable  size,  and  so  this  part  of  the  instrument  gap  remains. 

A proposal  to  use  an  instrument  of  the  above  type  in  a satellite  that  has  an 
equatorial  orbit  and  and  ranges  over  2 < L < 10  is  under  consideration  bv  NASA  as  part 
of  the  OPEN  (Origins  of  Plasma  in  the  Earth's  Neighborhood)  program.  The  other  part  of 
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the  program  calls  for  an  investigation  into  the  composition  of  the  interplanetary  plasma 
and  energetic  particle  population  by  constituents.  It  is  not  clear  that  the  two  missions 
are  compatible  on  a single  satellite,  since  the  orbit  required  for  the  interplanetary  phase 
would  (like  ISEE's  orbit)  be  unsuited  for  the  magnetospheric  phase  of  the  mission.  With  a 
well  chosen  orbit,  however,  the  OPEN  program  should  have  high  priority  as  a magneto- 
spheric  mission. 

A wide  variety  of  active  and  passive  plasma-injection  experiments  could  be 
performed  on  Space  Shuttle  flights  during  the  1980's.  The  low  altitude  of  the  Space 
Shuttle  orbit  will  necessitate  use  of  shaped  charges  or  Shuttle-launched  satellites  to 
inject  plasma  all  the  way  to  the  equator,  but  it  should  be  possible  to  inject  many  tens  of 
kilograms  of  ions  into  that  portion  of  the  magnetosphere  that  lies  outside  the 
plasmasphere,  into  the  dayside  cleft,  or  into  other  interesting  regions.  A general  picture 
of  what  can  be  accomplished  in  this  context  has  already  been  discussed  in  Sections  III  and 
V. 

As  a first  step  in  this  direction,  the  AMPTE  (Active  Magnetospheric  Particle 
Tracer  Experiment)  program  calls  for  injection  and  tracing  of  plasma  with  two 
conventional  satellites:  one  for  injecting  the  pies'-  . and  one  for  tracing  it.  Injection 
would  occur  not  only  in  the  magnetospheric  tail  but  also  in  the  solar  wind  near  the  nose 
of  the  magnetosphere,  using  lithium  and  perhaps  europium.  The  tracer  satellite  is 
proposed  to  have  an  equatorial  orbit  with  apogee  reaching  at  least  to  L = 8.  It  would 
carry  a set  of  composition  detectors  that  fill  the  instrument  gap  mentioned  above  and  in 
Section  I.  The  tracer  satellite  is  therefore  similar  to  the  OPEN  satellite;  either  could  be 
used  to  study  the  injected  ions. 

There  is  an  obvious  reason  for  using  tracer  ions  such  as  lithium  and  europium 
which  do  not  occur  naturally  in  the  magnetosphere:  they  are  analogous  to  the  red  dye 


injected  into  water  in  a hydrodynamic-flow  experiment.  As  long  as  one  is  interested  only 
in  studying  convective  transport,  there  is  no  difficulty.  However,  once  these  ions  acquire 
enough  energy  that  diffusion  outweighs  convection,  they  may  not  be  transported  as 
rapidly  as  protons  of  the  same  E/A  and  will  not  directly  reveal  what  happens  to  protons. 
Of  course,  such  heavy  ions  will  be  extremely  important  as  indirect  probes  of  the 
mechanism  responsible  for  radial  diffusion  in  general  (and  of  protons  in  particular)  when 
taken  in  conjunction  with  the  theoretical  framework  outlined  in  Section  III.  Such  radial- 
diffusion  studies  are  quite  readily  done  with  naturally  occurring  ions  in  the  magneto- 
sphere. The  real  value  of  exotic  ion  tracers  is  to  study  particle  injection  into  the 
magnetosphere:  to  determine  what  percentage  of  the  ions  from  a given  location  in  the 
solar  wind  or  geomagnetic  trail  ultimately  enter  the  region  of  trapped  particles,  where 
they  enter  it,  what  might  be  the  relationship  of  the  plasma  sheet  to  the  radiation  belts, 
and  so  on. 

The  last  observational  topic  on  our  agenda  relates  to  ion-composition  studies  in 
other  planetary  magnetospheres,  notably  that  of  Jupiter.  To  study  the  Jovian  magneto- 
sphere from  a satellite  in  Jovian  orbit  is  a high-priority  goal  of  the  near  future;  it  has 
been  combined  with  an  entry  probe  into  Jupiter's  atmosphere  in  the  Jupiter  Orbiter  Probe 
(JOP)  program.  The  Orbiter  should  be  capable  of  studying  the  magnetosphere  for  a long 
period  of  time.  This  will  require  an  orbit  at  r » 15rj  if  the  vehicle  is  to  avoid  severe 
radiation  damage  to  its  instruments.  Fortunately,  the  satellite  Ganymede  has  an  orbital 
radius  of  about  15rj,  a coincidence  that  will  enable  the  Orbiter  to  conduct  detailed 
studies  of  this  largest  Jovian  satellite  and  its  interaction  with  Jupiter's  radiation  belts. 
Unfortunately,  direct  studies  of  Io  and  Europa  would  not  be  possible  unless  the  Orbiter 
were  sent  into  a much  lower  orbit  near  the  end  of  the  mission.  The  Probe  would  report 


on  the  chemical  and  physical  composition  of  the  atmosphere  in  much  the  same  way  as  did 
the  earlier  Pioneer-Venus  probes.  Both  the  Probe  and  the  Orbiter  would  carry  ion- 
composition  detectors,  but  those  on  the  Probe  may  be  limited  to  low  energies.  Aside 
from  furnishing  data  bearing  on  Jovian  magnetospheric  physics,  the  Probe  would  yield 
much  information  of  astrophysical  interest,  since  Jupiter  is  believed  to  have  a 
composition  similar  to  that  of  the  primordial  solar  nebula. 

The  challenge  of  the  future  in  magnetospheric  heavy-ion  physics  is  not  entirely 

observational.  The  presence  of  several  distinct  ionic  species  in  a plasma  or  radiation  belt 

serves  to  resurrect  many  important  theoretical  problems  that  had  been  solved  for 

hydrogenic  plasmas  only.  The  additional  constituents,  for  example,  tend  to  introduce 

additional  wave  modes.  This  phenomenon  has  been  studied  in  Vlasov  theory  for  ion- 
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acoustic  waves  ",  but  a wide  variety  of  plasma  wave  modes  and  their  instabilities 
remains  unexplored  in  this  context.  For  example,  drift  waves  (which  are  associated  with 
spatial  gradients  in  a plasma)  are  commonly  found  to  show  maximal  growth  at 
wavelengths  2ir/k1  comparable  to  the  mean  ionic  Larmor  radius  / ft  , or  more 
specifically  at  k ^ ~ 1,  as  does  the  collisionless  tearing  mode  (which  may  be 

relevant  for  the  geomagnetic  neutral  sheet).  What  might  be  the  effect  of  a second  ionic 
species,  having  a different  mean  value  of  v^/n  , on  such  instabilities?  Would  this 
create  additional  wave  modes,  or  would  it  hybridize  the  modes  that  already  existed  with 
one  ionic  component?  Is  the  peak  growth  rate  enhanced  or  diminished  by  the  inclusion  of 
heavy  ions?  The  answers  to  such  questions  require  a serious  reconsideration  of  drift 
waves  and  tearing  modes  in  the  multi-ion  context.  The  effects  of  heavy  ions  on 
electrostatic  cyclotron  waves  are  similarly  unknown  and  in  need  of  investigation,  as  are 
the  relationships  between  electromagnetic  cyclotron  waves  involving  magnetospheric 
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heavy  ions  (such  as  0+)  and  geomagnetic  pulsations  (such  as  Pc-2  and  Pc-3)  observed  on 
the  ground.  It  should  be  clear  from  these  questions  that  much  work  remains  to  be  done 
on  multi-ionic  plasma  theory,  and  that  the  results  of  such  work  will  be  directly  relevant 
to  observations  made  in  space. 

Thus,  we  conclude  our  present  review  on  an  optimistic  note.  Heavy  ions  play  an 
important  diagnostic  role  in  magnetospheric  physics,  but  they  also  take  part  actively  in 
the  dynamics  of  the  mag letosphere.  They  are  suitable  for  magnetospheric  injection  as 
tracers,  but  they  can  also  be  used  to  trigger  conspicuous  upheavals  of  the  magnetospheric 
plasma.  Heavy  ions  are  important  constituents  of  the  earth's  magnetosphere,  but  they 
may  be  even  more  important  in  Jupiter's.  The  future  of  magnetospheric  heavy-ion 
research  (observational,  experimental,  and  theoretical)  is  indeed  a bright  one. 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 


Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,*  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine 


Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting. superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 


Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  oi  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor 
rosion  and  fatigue -induced  fractures  in  structural  metals. 


Space  Sc i. nee.  Laboratory:  Atmospheric  end  iono.ph.ric  phy.ic.,  radia- 
tion [roiti  the  atmo.ph.re,  d.n.ity  and  eompo.ition  of  the  atmo.pher.,  aurora, 
.nd  air.il>.  ma.neto.pheric  phv.ica,  cosmic  raya.  ten. ration  and  propafatton 
•f  plasma  —are.  in  the  ma. net o.ph* re . .alar  phv.ic.,  studies  of  solar  mafnatic 
A. Id.  ep.ee  sMtoMtnt,  a - rsy  astronomy,  the  effects  of  nuclear  euploaions. 
m stn.fi r storms,  and  snlsr  srtirtty  on  the  earth's  atmosphere.  I— ao sphere,  and 
mataeiassher.  the  .fferte  of  apei  rat.  .lectr  imapnatir , and  parti  relate  rndta- 


